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ABSTRACT 


Results of theoretical and experimental investigations of bellows 
typical of those found in Space Shuttle external tanks are presented. New 
correlation parameters are identified which generalize the alternating 
stress calculations presented in an earlier SwRI study titled "Bellows 
Flow-Induced Vibrations and Pressure Loss. 1 ' Alternating stress amplitudes 
and mean stress levels fonr. the basis of a fatigue analysis incorporating 
seven-ordinate charts for 347 S.S., Alloy 21-6-9, and Inco 718. A crack 
propagation model is also presented. Computer programs for computing 
bellows fatigue life and Two Phase flow and material hardness topics are 
contained in the report. 
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I. INTRODUCTION 


1.1 Overview 

This report describes all work performed by Southwest Research Institute 
under Contract NAS8-31994, "Research Study of Flow Induced Vibrations." This 
study was performed for George C. Marshall Space Flight Center of the National 
Aeronautics and Space Administration, and it was administered by the Structures 
and Propulsion Laboratory, with Mr. R. H. Veitch serving as Technical Manager. 

The general objective of this study was to evaluate bellows related 
theoretical assumptions either by analytical and/or experimental investigations. 
Emphasis was placed on obtaining a better understanding of the fluid-elastic 
excitation mechanism and upon developing a refined fatigue prediction methodology. 
The foundation of the current study is found in earlier research work performed 
by the Institute which is reported in a document titled "Bellows Flow-Induced 
Vibrations and Pressure Loss," by C. R. Gerlach, et al.^ 

Summary of Results 

A number of significant findings have been made throughout this report; 
these are summarized below. 

(a) Definition of Cp* Parameter - A stress correlation parameter has 
been defined which generalizes the existing bellows data contained 
in Reference 1. Previous data were characterized by a number of 
parameters such as the specific spring rate, fluid state, geometric 
factors and a vortex force coefficient. All of these factors are 
accounted for in the Cp* correlation and its usage. 

(b) Damping Model - As an alternate method of predicting stress ampli- 
tudes, an empirical damping model was developed. 

(c) Fatigue Prediction - A stress analysis has been coupled with the 
flow- induced vibration analysis in order to determine, with reason- 
able accuracy, the bellows fatigue life under varying environmental 
factors. 

(d) Computer Program - A computer program has been developed to allow 
quick computation of the bellows mode frequencies, lock-in ranges, 
stress Indicator, and stress level. 
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(e) Acoustic Resonance - The acoustic resonances as identified by 
analysis have been verified by limited experimental investigation. 

(f) Special Problems - During the course of the contract, several urgent 
and special bellows related problems were addressed at NASA's 
request. The solution of these problems are Included in this report. 

1.2 Organization o f Study 

The bellows study has been broken into two separate methods of approach 
as indicated in the block diagram shown in Figure 1. The end objective of both 
methods is to predict the fatigue life of U-shaped bellows made of an arbitrary 
material, and in both cases, the alternating stress component is generated by 
flow induced vibrations. Method I incorporates the stress indicator concept, 
while Method II Incorporates actual stress predictions which may be incorporated 
with 7-ordinate fatigue curves to predict bellows life. Method I suffers from 
the lack of a fatigue data base which must be generated by failing numerous 
bellows while influenced by flow induced vibrations. Method II suffers from 
underdevelopment of a realistic stress-deflection model where the convolute 
deflections can be predicted given an arbitrary geometry and flow conditions. 
Method I has been streamlined and somewhat generalized with the development 
of an envelope parameter designated as Cp* which is then used to determine the 
stress indicator. Method II efforts were directed toward the development of a 
flow induced stress model. 

1.3 General Discussion of Study 

The main propulsion system of the Space Shuttle is configured with three 
engines, a complex array of liquid and gas flow lines, and two large external 
tanks (ET). An elementary schematic of the main propulsion system is shown 
in Figure 2. Bellows are contained throughout the flow network; however, the 
bellows of primary Interest are contained in the feed lines (LO 2 and LH^) a nd 
in the small recirculation lines. 

Earlier studies have shown that unshrouded shuttle application bellows 
(see Figure 3 for bellows nomenclature) will vibrate violently when the contained 
fluid is moving at a specific critical velocity. The oscillation is shown to 
occur at a reduced velocity (U/fC) of approximately 4.5. Vortex shedding from 
the individual convolutes was found to be the flow induced vibration mechanism. 
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N • NUMBER OF CONVOLUTIONS COUNTED 
C FROM THE OUTSIDE 

N ■ NUMBER OF PLYS 
P 

D_ ■ MEAN BELLOWS DIAMETER 
m 

t • WALL THICKNESS (THICKNESS PER 
PLY IF MULTI -PLY) 

\ - CONVOLUTE PITCH 

0 * CONVOLUTE WIDTH 

a • MEAN FORMING RADIUS 

h • MEAN DISC HEICHT 


FIGURE 3. BELLOWS NOMENCLATURE 
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Experimental data, obtained from the earlier studies, were parametrically 
correlated in terms of (1) the Strouhal number vconvolute width is the charac- 
teristic dimension), (2) the bellows modal frequencies which Included added fluid 
mass terms, and (3) a stress Indicator which is proportional to the tiaxlmum 
dynamic stress. 

It has been shown that the stress indicator is a function of a vortex 
force coefficient, Cp, and a forced response dynamic amplification factor, Q. 
These experimentally derived factors are shown in Figures 4 and 5 and Table I. 

Finally, the observed fatigue life was related to the stress indicator 
as shown in Figure 6. The fatigue data were obtained for 321 S.S. only; 
although the general presentation could be expanded to include other materials 
if appropriate material factors could be included. Bellows for Space Shuttle 
applications are constructed of Inco 718 and steel alloy 21-6-9 Materials. 
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FIGURE 4. SUMMARY OF BELLOWS VORTEX FORCE COEFFICIENT EXPERIMENTAL DATA 






Curve Numbers Correspond To Certain 
Applications - See Table £ For These 
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FIGURE 5, DYNAMIC AMPLIFICATION FACTORS FOR VARIOUS BELLOWS APPLICATIONS 
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TABLE I. 

APPLICATIONS INFORMATION FOR USE WITH Q VALUE 
DATA IN FIGURE 5 


Specific Spring 

Number 

Internal Media 

Curve 

Rate (see Note 1) 

Plies 

(see Note 2) 

No. 

all ranges 

1 

low pressure gases 

i 

over 2000 lb/in^ 

1 

high pressure gases, light liquids 

i 

over 2000 

1 

water, dense liquids 

2 

under 2000 

1 

high pressure gases, light liquids 

2 

under 2000 

1 

water, dense liquids 

3 

over 3000 

2 

all 

3 

2000 - 3000 

2 

all pressure gases 

4 

under 2000 

2 

all pressure gases 

5 

2000 - 3000 

2 

all liquids 

5 

under 2000 

2 

all liquids 

6 

over 3000 

3 

all 

4 

2000 - 3000 

3 

all 

5 

under 2000 

3 

all pressure gases 

5 

under 2000 

! 

3 

all liquids 

6 


Use of Table - To use table, first calculate bellows specific spring rate, 
then look up application curve r Tiber corresponding to 
this specific spring rate, numbei. of plies, and internal 
media. 


Note It 


Note 2: 


The specific spring rate is here defined as 

K A N c 

^ « n _ A U 


or is the spring rate per convolute, per ply, per unit of 
diameter. 

Low pressure gases will be defined here as being those 
gases below 150 psia. Light liquids will be defined as 
having a density, relative to water, of less than 0.2. 




Stress Indicator - ( CfCppaQ/Np ) ( h/t ) x 10 psi 
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1.4 Review of Relevant Literature 

The following list of reviewed sources of bellows information is in- 
cluded to help direct the interested reader build a background knowledge which 
is needed for detailed evaluation of bellow related topics. 

(1) Kleppe, S. R. , "High Pressure Expansion Joint Studies" 

ASME Petroleum Mechanical Engineering Conference, New 
Orleans, Sept. 25-28, 1955, Paper No. 55-PET-10. 

A semi-torus expansion joint was extensively strain-gaged and 
hydrostatically tested. Test results compared favorably with 
R. A. Clark's theory as presented in "On the Theory of Thin 
Elastic Toroidal Shells," Journal of Mathematics and Physics, 

Vol . 29, 1950, pp. 146-178. 

(2) Turner, C. E. , and Ford, H. , "Stress and Deflection Studies of 
Pipeline Expansion Bellows," Proceedings of the Institute of 
Mechanical Engineering, pp 596-552, Vol. 171, No. 15, 1957. 

This paper presents a theoretical solution together with an 
experimental study of axial compression of certain bellows 
mainly of the corrugated-pipe type, used in the pressureless 
state. The total strain energy is written in terms of the 
circumferential stress and the axial loading moment. A 
Rayleigh - Ritz method is used to solve for a minimum strain 
energy condition. Ultimately the surface stresses are ana- 
lytically determined. The paper contains a short literature 
review covering the period from 1916 to 1953. 

(3) Feely, F. J., Jr., and Goryl, W. M. , "Stress Studies On Piping 
Expansion Bellows," Journal of Applied Mechanics, Paper No. 
44-APM-22 . 

In this paper a formula has been derived to show the total 
stress induced in the material as a result of the combined 
effects of pressure and movement. The validity of the approxi 
mations used in the formula have been verified by laboratory 
strain measurements. The paper deals primarily with flat disc 
type bellows. 
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(4) Samans, Walter, "Endurance Testing of Expansion Joints," ASME 
Paper No. 54-A-103. 

This paper presents the results of testing 19 bellows of various 
types to their endurance limit. The types include (1) welded 
roots, (2) hydraulically formed, and (3) welded disk. The 
bellows material consisted of stainless steel types 304, 321, 
and 347. A typical stress-distribution diagram for a 12-inch 
diameter hydraulically formed bellows is presented (case of axial 
extension and compression, and internal pressure). Strain mea- 
surements were taken with SP.-4 strain gages. The maximum stress 
range for both radial and circumferential stresses occurs near 
the root of the corrugation. 

(5) Haringy, J. A., "Instability of 3ellows Subjected to Internal 
Pressure," Philips Res. Report 7, 189-196, 1952. 

Jellows may become unstable when loaded by internal pressure. 

The critical value of this pressure is governed by the rigidity 
of the bellows with respect to bending. Critical pressures have 
been analytically determined for rectangularly shaped corruga- 
tions and these critical pressures may be considered to agree 
with those obtained experimentally for U-shaped bellows when 
considering the approximations introduced and the variation of 
wall thickness. 

(6) Laupa, A., and Weil, N. A., "Analysis of U-Shaped Expansion 
Joints," Journal of Applied Mechanics, Transactions of the ASME, 
pp 115-123, March 1962. 

An elastic analyses of U-shaped expansion joints under axial 
loads and internal or external pressure is presented. The 
analysis employs the energy method for the toroidal sections, 
and the theory of symmetrical bending of circular plates 
augmented by thick walled cylinder analysis for the annular 
plate connecting the two toroidal sections. 

(7) Sack, L., "Avoiding Fluid-Line Failure in Bellows and Convoluted 
Tubing," Machine Design, May 27, 1971. 

Flexline response frequencies are modeled as a lumped parameter 
system where the characteristic frequency is determined by 
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knowledge of the convolute effective mass and the effective 
fluid compressibility. Bellows longitudinal natural frequencies 
are modeled as a spring-mass analog where a dimensionless fre- 
quency parameter is utilized for evaluating all the longitudinal 
modal frequencies. An attempt has been made to define the 
maximum alternating stress. 

(8) Baylac, G. , et al., "Calculation of Acoustical Resonances in 
Irregular Cavities with Application to Noise-Induced Stress in 
Expansion Joints," ASME Paper No. 75-DET-64. 

An analytical and experimental study of the acoustic behavior 
of seven and nine corrugation expansion joints (bellows) used 
in a nuclear reactor is presented. Resonant frequencies obtained 
from a computer program using a matrix method are given. Experi- 
mental test results on seven corrugation expansion joints are in 
good agreement with the computations. It is concluded that the 
calculation of acoustic frequencies of expansion joints with in- 
ternal sleeves can be utilized to avoid the coincidence of these 
frequencies with those of a mechanical or flow-induced noise 
nature and thus reduce the loads on expansion joint corrugations. 

(9) T. M. McCrary, "Evaluation of Inconel 718 Bellows Material," 
SD73-SA-0014, Rockwell International Space Division, Mar. 1973. 

Life cycle testing was performed on 10" diameter bellows with 
nominal 3/8-inch high convolutions (.008-inch thick. Inconel 718). 
Testing was similar to that conducted for Boeing Company by 
Strazar. Metallurgical and fatigue properties were evaluated. 

This report does present a source of fatigue data as a function 
of bending stress (bellows), and percent of tensile ultimate 
strength (specimens only) . 

(10) "Effect of Surface Irregularities on Bellows Fatigue Life," 

R7250 Rocketdyne, NASA Contract NAS8-19541. 

The report presents the results of a brief test program aimed at 
generating data on bending life of notched CRES sheet specimens. 
Emphasis of the study is directed toward the quantitative nralu- 
ation of bellows' defects, particularly those resulting from 
accidental damage. An empirically derived procedure for evalu- 
ating bellows' surface irregularities and determining service 
life is presented. 


4 
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II. GENERALIZED CORRELATION PARAMETER: C F * 

11. 1 Introduction 

Through the efforts of Gerlach, et al.^'J and Sack^®) is has been 
well established that a series of lumped spring-mass elements can represent 
a free bellows and the modal frequencies can be computed with a high degree 
of accuracy. The work of Gerlach went on to show that the flow excitation 
mechanism is a vortex shedding phenomena that occurs in the entrance region 
of a convoluted bellows. When the vortex shedding frequency is near a 
bellows longitudinal structural frequency, the vortex shedding frequency 
will "lock-on" and the structure will vibrate at an amplitude dependent 
upon the amount of fluid and structural damping present. 

Ultimately, the most fundamental question is how to determine the 
amplitude of convolute displacement and hence the resultant maximum alter- 
nating stress amplitude. Two stress prediction models wi'l be addressed 
in this section. 

11. 2 Cp* Correlation Parameter 

Reference 1 contains the derivation and application of a stress indi- 
cator concept. It must be emphasized that the original form of the stress 
indicator was merely a bench mark showing relative stress intensities as a 
function of fluid and geometric parameters. Its purpose was to guide a de- 
signer when obtaining fatigue predictions. The stress indicator concept is 
a valid method for predicting fatigue life so long as a substantial data base 
is developed; unfortunately, a large data base does not exist. 

Before describing the Cp* ("C sub F Star") model, the original stress 
indicator model is reviewed. It has been shown that the maximum convolute 
stress due to flow induced vibration is 

<?ait - K ^ (h/t) 2 1/2 p f V crit 2 (II. 1) 

The K term contains factors of proportionality relating to geometric 
constraints and this factor was extracted from Equation (II. 1) to produce a 
single simple expression for stress which contains only readily known bellows 
dimensional data, parameters, and flow variables. Therefore, the indicator 
is given as 
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S.I. - (h/t) 2 1/2 p f V crit 2 (II. 2) 

Table II compares the calculated stress indicator and measured stress 
on the crown of the second convolute (see Appendix B for a description of the 
experimental techniques). Several items are worth noting in this table. The 
K factor ranges from 0.585 to 3.61 for the limited test conducted and there 
is a downward trend in the K factor as the mode number increases. This shows 
that the stress indicator may or may not be a conservative estimator of stress 
levels, and the K factor is not constant as assumed in Reference 1. 

TABLE II. MEASURED CONVOLUTE RADIAL STRESS AND CALCULATED 
STRESS INDICATOR COMPARISON 


Bellows* 
Ident . 

r ,-, — 

Mode 

No. 

Measured 
Radial Stress 
KSI (peak) 

Stress 

Indicator 

KSI 

Measured 

Calculated 

4 

— 

i 

2 

2.21 

1.325 

4 

2 

8.02 

8.24 

.97 

4 

3 

8.94 

11.48 

.775 

6 

1 

.765 

.93 

.82 

6 

2 

3.67 

3.70 

.99 

6 

3 

4.59 

7.83 

.585 

15 

1 

2.82 

.78 

3.615 

1' 

2 

7.84 

3.18 

2.46 

15 

3 

8.51 

6.76 

1.225 

E 

1 

1 

4.57 

1 

2.28 

2.00 

E 

2 

8.41 

7.91 

1.05 

*Dimensional Data 

is contained in Appendix C. 



The stress indicator contains two terms, Cp and Q, that are 
dependent upon factors of damping, internal pressure, convolute geometry, 
and the flow media. Values for Cp are obtained from Figure 4 while values 
for Q are obtained from Figure 5 and Table I. The data contained in these 
sources have been correlated in the form of one universal stress function 
curve as discussed below. 
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All data contained In Reference 1 has been evaluated in terms of a 
correlation parameter defined as 

C F * - C f Q(N/N c ) 

Figures 7 through 9 show plots of the force coefficient parameter for 
representative samplings of the total data base. The effect of changes in X 
on the force coefficient parameter Cp* is illustrated in Figure 7. Here, a 
single bellows was tested at various pitch values, X, and the peak response 
of the first longitudinal mode (N-l) was noted. It is noted that spring rate 
is affected somewhat by changes in X. 

The reduced data shown in Figure 8 clearly illustrates the effect of 
vortex reinforcement and vortex retardation on the flow induced response of 
the bellows. 

A vortex reinforcement occurs when the vortex shedding from an upstream 
convolute arrives at the adjacent downstream convolute at the right moment 
to aid in the formation of the vortex forming at that adjacent convolute. 
Vortex retardation has the opposite effect. The vortex shed from an up- 
stream convolute arrives at the adjacent downstream convolute at the right 
moment to detract from the formation at that location. As we will soon 
discuss, it is our present concept that vortex reinforcement is most pre- 
valent and effective in the higher longitudinal modes. (Figure 6 from the 
final report "Bellows Flow- Induced Vibrations and Pressure Loss" clearly 
shows a visualization of vortex reinforcement for a higher longitudinal 
mode.) In the first two or three modes of a bellows, vortex reinforcement 
and vortex cancellation both come into play, as illustrated by Figure 8. 
However, for the intermediate modes, the vortex retardation phenomena is 
prevalent . 

Figure 8 presents a plot of Cp* versus the mode number N for four 
test bellows that have constant values of the parameter (h/t) but have 
h values ranging from 0.2 to 0.5. Since spring rate is proportional to 
(h/t), this family had similar modal frequencies, so that the effect of 
convolute height h should be revealed. Also, however, each of the four 
bellows was tested for three or four values of X achieved by stretching. 

Note ''hat there is a spread of the combined Cp* values for these four 
bellows for each mode number of N value. This spread is caused by a 


C F 0 (N/N 
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A/a 

1.5 2.0 2.5 



FIGURE 7. VORTEX FORCE COEFFICIENT Cj* VS. PITCH FOR 
THE FIRST MODE OF BELLOWS 105 
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FIGURE 8. VORTEX FORCE COEFFICIENT C F * AS A FUNCTION OF 
MODE NUMBER FOR BELLOWS WITH CONSTANT (h/t) 
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combination of two factors. First, it represents the Influence of the effect 
of changing X as illustrated previously in Figure 7, and, secondly, it re- 
flects the normal variation expected in flow-induced vibration experiments 
of bellows where slight changes in alignment, clamping of the ducting, etc. 
cause changes in the peak response point. 

From Figure 8, we have concluded the following: 

(a) Other than for the No. 1 specimen, which had h*0.2 or a 
very short convolute, the effect of h was not apparent 
between the bellows. Specimen No. 1 had lower Cp* values 
than the other bellows, probably because short convolutes 
do not couple so well as taller convolutes. After all, the 
limiting case is h • 0 which represents a straight pipe 
which has no response of the type under consideration. 

(b) The vertical spread of Cp* for each mode is primarily caused 
by vortex reinforcement or vortex cancellation. 

(c) The pronounced minimum of Cp* is a result of an optimum 
vortex cancellation effect for this mode number range. 

(d) The rapid rise of Cp* for the higher longitudinal modes is 
a result of a predominance of vortex reinforcement for 
these modes. 

(e) Many of the higher modes simply never appear because other 
modes close to them predominate and prevent their occurrence. 

Figure 9 presents Cp* as a function of mode number N for three bellows 
having similar convolute geometry but different numbers of convolutes. The 
bellows No. 19 illustrates yet another phenomena. Note that tne Cp* values 
for this bellows are quite low for the first two longitudinal modes. Aslo 
note the strong presence of the first cocking mode plotted for N * 1.5. 

For this bellows the cocking mode was stronger than normal so it suppressed 
the first and second longitudinal modes causing their Cp* values to be ab- 
normally low. 
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The primary intent of tne Cp* relation is to mathematically collapse 
all of the experimentally generated Q surfaces into one relationship that 
applies to all ranges of the bellows operational parameters; hence, the 
stress indicator is computed 

C * N 

S.I. * -| ~ (h/t) 2 (1/2 p V crit 2 ) (II. 3) 

The parameter Cp* is obtained from Figure 10 which is a somewhat 
conservative curve that envelops all previously generated experimental 
bellows data. This curve contains all inherent information relating to 
Cp and Q . 

II J Summary of Design Analysis Procedure 

The procedure for analyzing a given bellows design to assure freedom 
from flow-induced vibration failure consists of several distinct steps which 
are listed below. 

Step 1. Calculate the natural frequencies for all modes of the 
bellows . 

Step 2. Determine the lock-in or critical velocity range for each 
possible mode of vibration. 

Step 3. Calculate the Stress-Indicator for each mode at the critical 
velocity. 

Step 4. Determine the potential for failure of the bellows using the 
Stress-Indie, tor versus Cycles-to-Failure curve. 

Pages 23 and 24 r resent a detailed step-by-step procedure that may be 
used for hand calculations. A more sophisticated calculation procedure is 
contained in a computer program (see Appendix A). 
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FIGURE 10* ENVELOPE CURVE FOR C F * CORRELATION 





STEP A 
STEP B 


STEP C 
STEP D 


STEP E 


TABLE III. 

SUMMARY OF FREQUENCY AND STRESS LEVEL CALCULATIONS 


23 


Consider the bellows structure representable by a lumped 
mass-spring mechanical model. 


Calculate the elemental spring rate value K from the 
expression 


K - 2 N c K a 

where K A is the overall spring rate determined from a force- 
deflection test or from the following expression: 

K A * Dm E ^ (t/h ) 3 


Calculate the elemental metal mass 


M m - ffProt N p D m [TTa + (h-2a)] 


Calculate the fluid added mass Mf, for the first few 
longitudinal modes and for the higher longitudinal modes as 


and 


M f * tt/ 2 p f D m h (2a-tN p ) 


M f 


* P m Pf h ' 
36 


^ First few N values 


^ Higher N values 


Calculate the reference frequency f Q from the expression 


f Q - 1/2 tt y/k/m 


where m * 1 % + mf 


TABLE III (CONTD) 

SUMMARY OF FREQUENCY AND STRESS LEVEL CALCULATIONS 
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STEP F 


STEP G 


STEP H 


STEP I 


Calculate the dimensionless frequencies and then multiply the 
dimensionless frequencies by the reference frequencies to 
obtain the true mode frequencies 

_ /_ I, , /ir(2N c -i)\| 1 Dimensionless frequency 

i ”V L C0S V 2»c ) \ J for the i-th mode 

1 - 1,2,3, ... 2N C -1 

fi * f 0 \ True frequency for the 

J i-th mode 

Alternately, the dimensionless frequency factors may be 
obtained from Table I» Appendix A. 

Calculate the first convolute bending mode from the 
expression 


f b * l/2ir \J 8k/m 

where m = + m£ 

and * tt D m p£ h 3 /36 

Calculate stress indicator from the following expression: 

s.i. - C F * (^) (h/t) 2 a/2 p V crlt 2 ) 

The parameter Cp* is obtained from the curve presented in 
Figure 10. 


Calculate bellows expected life from the data presented in 
Figure 6, which is a plot of stress indicator versus cycles 
to failure. If the fatigue life is greater than 10 s cycles, 
then the data are conservative for materials classified as 
Inco 718 and alloy 21-6-9. 

If the calculated number of cycles is less than 10 5 , then 
the expected life cf alloy 21-6-9 will be less than that 
indicated for SS-321 or its equivalent SS-347. 
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III. STRESS LEVELS 


111.1 Introduction 

While section II presented a method for calculating vibration 
frequencies and stress-like quantities that may be used with the appropriate 
analysis to predict fatigue life, this section will explore various properties 
of actual stress levels experienced during the flow induced vibration process. 
As of this writing, an exact method has not been developed to calculate actual 
stresses; however, several important aspects of the problem are presented 
along with a reasonable stress :alculation procedure. 

111. 2 Stress Envelope 

Test data, shown in Table IV, has been reauced in terms of non- 
dimensional stress and velocity ratios for each longitudinal mode of vibra- 
tion. The velocity ratio is formed by dividing the critical velocity of a 
particular mode by the first mode velocity and the stress ratio is formed in 
a similar fashion. The correlation in Figure 11 shows that similar families 
of curves are developed. The data may be further collapsed by referencing 
the curves to a particular damping ratio. For the present case, an average 
damping ratio of .00635 served as the reference damping value. Figure 12 
shows the results of the damping normalization. From the limited data pre- 
sented, the second and third mode stress may be calculated by the following 
empirical equation, 



alt 




(III.l) 


where F£ * 2.75 
F 3 * 3.05 

Equation III.l was developed from data obtained from a series of 
3", 321 S.S. bellows with a constant convolute height. The material thick- 
ness, number of plys and number of convolutes were allowed to vary and the 
measured spring rates were significantly different. The alternating stress 
component referred to is the convolute radial stress. Radial stresses were 
calculated from biaxial strain data (radial and circumferential) as described 
below. 
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Each bellows was strain gaged (see Figure 3, Appendix B) on the second 
and middle convolute in the radial and circumferential directions which are 
the assumed principal directions. Principal stresses are calculated from the 
measured principle strains. 

Or - j (e R + p e c ) (III. 2) 

l - V 


a 

c 



(£ c + V £ R> 


where 

a R * radial stress, psi 

o c * circumferential stress, psi 

E * modulus of elasticity, psi 

y = Poisson ratio 

Cr = radial strain, nucroinches 

e c ■ circumferential strain, microinches 


(III. 3) 


III. 3 Two-Ply Bellows 

Multi-ply bellows flow-induced strain characteristics are significantly 
different than those of single-ply bellows. Figure 13 shows the flow-induced 
strain for a 3" single-ply bellows and a 3" two-ply bellows. In each case, the 
first mode has been flow excited. Note that the alternating strain level 
for the single-ply bellows is independent of internal pressure, while the 
strain magnitude and lock-in range for the two-ply bellows is strongly dependent 
upon internal pressure. For the particular bellows exhibited, it was found that 
the alternating strain component varies inversely and as a linear function of 
pressure (see Figure 14). 


The most plausible explanation of this phenomena is that Coulomb friction 
damping is experienced between the plys of the bellows. The Coulomb friction 
force is directly proportional to the normal force acting in a manner to 
compress the plys together. To bear out this fact, a two-ply bellows was 
impulsed into vibration and then allowed to decay. The decay traces are 
shown in Figure 15 where it is obvious that the damping is a function of the 
internal pressure which is th mechanism generating the normal force on the con- 
volute sidewalls. 



VL-Strain, P-P 
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The results of these pressure tests suggest that multi-ply bellows 
vibrate with a lesser magnitude when they are internally pressurized; thus, 
when single-ply stress calculations are performed on a multi-ply bellows ex- 
hibiting the same damping ratio at zero gage internal pressure, the calculated 
alternating stress component will be conservative (higher stress) for the 
internal pressurization case. These tests suggest that it may oe practical 
to Include damping material between plys as an alternative to including flow 
liners. 

III. 4 Convolute Mean Stress 

Typically, alternating stres es which are generated by flow induced 
vbrations are superimposed upon a mean stress which results from internal 
static pressure and/or bellows axial extension or compression preload forces 
By observing a typical seven-ordinate fatigue chart (for example, see 
Figure 23), it is noted that fatigue life is decreased with increasing mean 
stress magnitude. For example, a bellows that is operated at high static 
pressures would fail sooner than one operated at lower pressure even if the 
alternating stress component were equal for both cases. The derivation and 
use of the seven ordinate curves will be discussed in Section IV; however, 
the important issue is that the seven ordinate charts allow for mean stress 
contribution which is not present in cycle-to-failure (S-N) curves. 

III. 4.1 Internal Pressure Stress 

Figure 16 presents the strain data obtained on Bellows No. SwRI-E 
during an internal pressurization test (ends of the bellows were clamped). 

The maximum principal stress was calculated for Convolute No. 2 and No. 7 
and these stress values compared to the following equation taken from 
Reference 4. 

a p - P/2 (h/t) 2 (III. 4) 

Table V summarizes the results which are evaluated at a pressure 
of 30 psig. Note, o p is a compressive stress on the convolute crown. 

The table compares the compressive stress, o p , to the measured radial 
stress, a max . 
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TABLE V. INTERNAL PRESSURE STRESS AT 30 FSI 


Convolute No. 

?max (KSI) 

Op (KSI) 

% Error 

2 

-24.9 

-21.1 

-18 

7 

-33.0 

-21.1 

-56 


It is noted that Equation (III. 4) under-predicts the radial stress 
(maximum principal) by as much as 56%. It is also noted that the radial 
stress in the center region of the bellows is higher. Most likely this 
higher center stress is caused by a "ballooning" effect in the mid-span of 
the bellows. The conservative approach when considering multi-ply bellows 
is to assume that the plys are not in complete contact; thus, the effective 
thickness is less than Np *t. Due to the limited data obtained with re- 
spect to ply-coupling effects, it is recommended that the calculated single 
ply stress be used when multiple plys are incorporated in a design. 

III. 4.2 Compression Preload Stress 

The same 6" bellows that was used for pressure tests was subjected 
to compression loading test. This is accomplished by placing calibrated 
weights on the open end edge of a free bellows which is placed in an upright 
position on a hard surface. This procedure is used to obtain the bellows 
spring constant K^; however, in this test the strain gage readings are 
also recorded. Figure 17 shows the strain data obtained versus compression 
loads. By noting that 



= change in microstrain per unit change, in live length 
(pc/ in) 

■ bellows spring constant (lb/in) (slope of deflect ■* un-load 
curve) 

* change in microstrain per unit change of load (ye/lb) 
(slope of strain-load curve) , 


where 


d ye 
di 


d ye 
d F r 




34.95 ye/lb ra = 36.66 pe/lb 



S3X°A a?PT-ia 


S3X°A aSptag 


FIGURE 17. STRAIN DATA FOR AXIAL COMPRESSION LOAD - 6" BELLOWS NO. 
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it is possible to determine the convolute strain-load characteristic. The 
deflection-load curve is pres' .id below (Figure 18) Trom which the bellows 
spring rate, K^, can be detenu id. 



FIGURE 18. LOAD-DEFLECTION CURVE - BELLOWS NO. E 

Axial compression stresses as obtained experimentally have been 
compared to the following equation: 

a c - (III. 6) 

h N c 

where 

O c « stress due to compression or extension load (a c > 0 for 
compression load, O c < 0 for extension load), psi 
A “ deflection of live length, inch. 

Table VI has been prepared to compare experimental results with Equation 
(III. 6) for a preload of 20 lbs. 
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TABLE VI. PRELOAD 

STRESS AT 20 POUNDS 


Convolute No. 

<W «si) 

O c (KSI) 

% Error 

2 

22.96 

22.09 

-3.9 

7 

24.21 

22.09 

-9.6 


It is observed that Equation (III. 6) gives reasonable accuracy and 
it provides a for relating relative convolute motion to convolute 

radial stress level. Equation (III. 6) can be used in a dynamic situation; 
however, it must be emphasized that the deflection value used is relative 
to adjacent convolutes. 

Equation (III. 6) is easily modified to incorporate preload rather 
than deflection if the bellows spring constant is known. 


thus , 


A = F C /K A 


Et F c 
h2 N c K A 


(III. 7) 


III. 4. 3 Compression Preload With Internal Pressure 

A schematic illustrating the nature of the radial fiber strains in 
the region of the bellows root and crown is shown in Figure 19. The strains 
are the result of bending moments generated in root and crown. For ana- 
lytical considerations, the bellows is envisioned to be restrained by the 
external piping for the case of internal pressurization. It is immediately 
obvious from Figure 19 that while it may be possible to reduce the crown 
radius stress state by simultaneous compression loading and internal pres- 
surization, the root stresses are intensified by the combination loading. 
Therefore, the root stress may be estimated as follows: 

o cp - 0 p + a c (III. 8) 

where a cp ■ combined stress due to pressure and compression load. By sub- 
stitution of Equation (III. 7) and (III. 4) into (III. 8). 
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). 

i. 

f 

i 



T 

T * Tension (+) C * Compression (-) 
(b) Compression Preload 


C 



T 

(c) Internal Pressure 
(Ends Restrained) 


FIGURE 19. NATURE OF FIBER STRAINS 
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■'cp 


P/2 (h/t) 2 + 


Et F c 
h 2 N c K A 


(III. 9) 


III. 5 Convolute Alternating Stress and Displacement 


A series of three-inch diameter bellows were flow tested to vali- 
date several assumptions made in earlier studies (Reference 1) . The vibra- 
tory peak stress in the bellows convolute was assumed to be given by 


a 


P 


c s Et X 


(III. 10) 


where C s is a geometric stress factor and the other terms are as defined 
earlier. The Reference 1 work utilized a single point strain gage to infer 
displacement and stress which is difficult under the best of test coriditions. 
In the present study, stress was measured via a biaxial gage arrangement and 
convoluted displacement was obtained independently via a displacement probe 
(see Appendix B for details) . 

By assuming a mode shape over the first quarter wavelength of the 

form 


X - X q /2 [(N/£)y -i- sin (N7ry/Jc)j (III. 11) 

where X denotes the axial absolute displacement of a given point along the 
bellows defined by the axial position coordinate y, we may determine the 
relative displacement by differentiating Equation (III. 11) with respect to y. 
Thus, 

A<S - X 0 /2 [(N/i) + Nir/fc cos (Niry/i.)] (III. 12) 

The above method was used to convert absolute displacement, 6, data 
into equivalent relative displacement, A6. 

A summary of the deflection and stress results are shown in Table 
VII for each test specimen at the first, second, and third modes and a 
summary of the damping characteristics is shown in Table VIII. 

Calculated alternating stress levels as determined by Equation 
(III. 10), have been correlated with actual measurements. Results shown in 
Figure 20 indicate that C s may be considered to equal unity. 
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TABLE VII. THREE-INCH BELLOWS DEFLECTION AND STRESS RESULTS 


Specimen 

No. 

Mode 

No. 

Vp, fps 

2nd Convolute 


A6, mills 

a Alt» ksi 

4 

1 

5.40 

4.0 

2.33 

2.93 


2 

10.80 

15.0 

3.21 

8.02 


3 

15.89 

22.6 

3.58 

8.94 

6 

1 

4.18 

0.8 

0.40 

0.765 


2 

8.35 

4.8 

1.65 

3.67 


3 

12.53 

7.8 

1.93 

4.59 

15 

1 

7.26 

4.0 

1.96 

2.82 


2 

14.52 

14.0. 

4.71 

7.84 


3 

21.79 

21.0 

3.34 

8.51 





























ksi 


42 



FIGURE 20. ALTERNATING STRESS VERSUS DEFLECTIONS 
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For all single ply bellows tested, the alternating stress component 
was observed to be Insensitive to Internal pressure variation. Tests were 
conducted at pressures of 0, 10, 20 psig over the first three modes of exci- 
tat Ion . 


III. 6 Recommended Stress Prediction Equation 
III. 6.1 Alternating Stress 

Two methods are available for calculating the alternating stress 
component. The Stress Indicator may be used as a predictor of actual stress 
for single ply bellows by incorporating a factor 2 into the S.I. equation 
over the first three modes of vibration, or 


°alt " 2 Cf * N c (h/t) 2 (1/2 p V crit 2 ) (III. 13) 

N N p 


The second method is merely a refinement of the above method. The 
stress envelope factor may be applied to the value of the S.I. calculated 
for mode 1, or 


3 alt 


N 


S.E. 


. 00635^ „ 

~T~) f n 


(III. 14) 


where F 2 * 2.75 

F 3 - 3.05 

The second method requires more detailed knowledge of the bellows; 
however, it provides a means to infer the effects of combined fluid and 
structure damping. 


III. 6.2 Convolute Mean Stress 

Significant errors have been observed in the measured and calculated 
stress values that relate to internal pressure while axial extension of com- 
pression preloads may be more accurately modeled. Therefore, the recommended 
mean stress model is 

0 mean * * d>/t) 2 + |4| (III. 15) 

where A is the total live length extension or compression displacement. 

For the multi-ply case, it has been assumed that t' e plys do not fully couple 
hence, the calculations for the single ply are applied to multi-ply designs. 
Due to the inaccuracy of the simple pressure-stress model, the .5 factor has 
been deleted. 
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III. 6. 3 Combined Stress 

The maximum stress developed Is composed of the two additive compo- 
nents, or 


^max “ ^alt + ^mean (III. 16) 

and by substitution of Equations (III. 13) and (III. 15), the proposed com- 
bined elastic stress model is 


°max 



•!- (h/t) 2 


P + 


2C f * N c 
N N p 



(III. 17) 


III. 7 Material Hardness Properties 

A 3” bellows with 13 convolutes was sectioned and prepared for micro- 
hardness testing. This is accomplished by cutting axial strips of approxi- 
mately 1/2" wide that contain several convolutes. Subsequently, these sec- 
tions are imbedded in an epoxy molding compound, then the compound and bellows 
specimens are ground until their surfaces exhibit a highly polished finish. 

The bellows specimen is placed into a Diamond Pyramid Hardness (DPH) testing 
apparatus where a specific sized diamond needed is allowed to penetrate the 
bellows surface. The driving weight used is 10 kg. By an appropriate mea- 
suring technique, the dimensions of the penetration, rhomboid shaped, are 
measured and then converted into a DPH number. 

Figure 21 shows a general bellows section. Three measurements were 
taken at the approximate locations shown in the figure; therefore, 24 data 
points per convolute were obtained. The results are tabulated in Table IX. 

Upon careful review of the data, several observations are apparent 
which include the following: 

1. Global averaging of the convolute center region produced a lower 
average hardness than found in global averaged outer edges. 

2. Zonal averaged hardness numbers in the "inside diameter" region 
exhibits hardness close to slightly below the global average. 

3. The "outer diameter" region exhibits hardness numbers signifi- 
cantly larger than the global averages. 

4. The "straight wall” region exhibits hardness numbers signifi- 
cantly lower than the global averages. 
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FIGURE 21. CROSS SECTION OF BELLOWS SHOWING LOCATIONS 
OF HARDNESS MEASUREMENTS 


LOCATION 


TABLE IX. HARDNESS READINGS 

DIAMOND PYRAMID HARDNESS 



CONVOLUTIONS 1, 
OUTER SURFACE CENTER 

2, & 3 

INNER SURFACE 

h l 

224 

213 

225 

32 (inside diameter) 

2'<1 

235 

246 

°2 

232 

242 

242 

C 2 (straight wall) 

239 

217 

232 

d 2 

268 

272 

264 

e 2 (outer diameter) 

270 

268 

284 

h 

266 

268 

258 

g 2 (straight wall) 

239 

224 

231 

h 2 

241 

226 

241 

a e (inside diameter) 

235 

218 

235 


mean (y ) 

Std. Dev. (a) 


245.5 

16.35 


240.1 

21.5 




CONVOLUTIONS 6 , 7, & 8 

h 6 

236 

263 

a-j (inside diameter) 

236 

253 

b 7 

257 

268 

C 7 (straight wall) 

235 

250 

d 7 

281 

2/9 

e-j (outer diameter) 

273 

265 

f 7 

285 

279 

g 7 (straight wa 1 !) 

236 

232 

h 7 

281 

253 

ag (inside diameter 

265 

257 


b 8 


268 

250 

273 

mean 

(y) 

259.4 

259 

264.3 

Std. 

Dev. ( 0 ) 

20.3 

13.78 

19.15 
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TABLE IX. 

HARDNESS READINGS 

(Cont 'd) 


LOCATION 

DIAMOND PYRAMID 

HARDNESS 


CONVOLUTIONS 11 

, 12, & 13 


OUTER SURFACE 

CENTER 

INNER SURFACE 

h ll 

236 

230 

236 

a^2 (inside diameter) 

221 

247 

263 

D 12 

233 

239 

261 

c^2 (straight wall) 

247 

236 

243 

1 d 12 

275 

275 

273 I 

1 

; e^2 (outer diameter) 

313 

275 

290 j 

f 12 

300 

294 

285 1 

1 

j gi2 (straight wall) 

243 

236 

236 

| h 12 

267 

265 

267 ! 

t 

| a]j (inside diameter) 

261 

265 

255 

bi3 

275 

268 

257 

j - 

mean (y) 

261 

257.3 

260.54 

Std. Dev. (0) 

286 

20.73 

17.96 
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5. Outer diameter zonal averaged DPH numbers ranged from 274 to 2 93 
which corresponds to a Rockwell Hardness range of 26C to 29C. 

6. Inner diameter zonal averaged DPH numbers ranged from 229 to 26C 
which corresponds to Rockwell readings in the range of 96B to 24C. 

III. 8 Conclusions 

1. The outer diameter region exhibits a yield stress of approxi- 
mately 132,000 psi whereas the inner diameter region exhibits a 
yield stress of approximately 100,000 psi. These yield values 
are somewhat lower (3C7) than those reported in Reference (13); 
however, it is speculated that the hydroforming process work 
hardens to a lesser extent than the rolling process. 

2. Failures most often occur in the root or crown region; there- 
fore, in view of the hardness data, it can be concluded that 
failures are not the result of material weakness in the failure 
region. 
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IV. FATIGUE LIFE 


IV. 1 Crack Propagation Model 

A bellows fatigue life model was developed based on the assumption 
that crack propagation in the convolute wall is the failure mechanism. It 
was further assumed that the crack was initiated by a pre-existing surface or 
material flaw. The state of stress in the bellows wall was taken to be the 
sum of the mean stress due to internal fluid pressure plus a cyclic bending 
stress that is associated with convolute deflection in any given mode of 
excitation. The stress model used here is different than that used in Sec- 
tion III; however, the features of the crack modeling and general results are 
valid. 

(12 13 } 

The mean internal pressure stress is, ’ ' 

a p = p/2 (h/t) 2 (IV. 1) 


where 

p 

= internal pressure. 



h 

= root-co-crown height. 

and 


t 

■ bellows wall thickness 

/■"■N 

H* 

X 

rt 

X) 

M 


Superimposed onto this steady state stress is a cyclic, deflection- 
related bending stress that is caused by the flow induced vibration of the 
bellows convolutes at given excitation mode. The peak-to-peak amplitude of 
this cyclic stress component is given by,^ 2 >13) 


A a 


2(1.5) Et 
A 1 2 



(IV. 2) 


where E = Young's modulus 
A = convolute pitch 

A/N c * flow- induced convolute deflection 


The deflection per convolute is calculated from Reference 1, 


where 


P 

V 

A P 

D i» D 0 

KA 

CpQ 

*“m 


_A_ 

N c 


82 K a 


(CpQ) 


(IV. 3) 


fluid density 

critical flow velocity as a function of mode number 
tt/ 2 h (Di + D 0 ) 

bellows inside and outside diameter, respectively 
bellows spring rate 

force amplification factor from Figures 4 and 5 
bellows mode factor. 
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The bellows mode factor. 


C ffl , la of the form^) 


Cm 



(IV. 4) 


where N * mode number 

N c ■ number of bellows convolutes. 

Thus, Equations (IV. 1) through (IV. 4) defiue the mode-dependent 
state of stress in the bellows wall . This state of stress can be Illustrated 
schematically as shown in Figure 22. In this figure, tensile stresses are 
positive. Depending on the mode number and the magnitude of the mean stress, 
the minimum stress can be compressive, in which case the sign of the stress 
is negative. 

If a crack is initiated on the bellows surface, the rate at which 
it will propagate into the wall thickness is governed by 


da „ ( YAa/"a ^ m 
dn * C \~I-R / 


where a 
n 

A0 

C,m 


crack length 

number of imposed stress cycles 
cyclic stress range. Equation (2) 

curve fit coefficients that describe the experimental 
crack growth rate as a function of stress intensity 
factor, which is the expression within the brackets in 
Equation (IV. 5). These parameters are dependent upon the 
bellows material. 


The factor, R, accounts for the mean stress effect, and jt is defined as 


Op-Ao/2 

cjp+Ao/2 


(IV. 6) 


It is worth noting that when Op is equal to zero (no pressure 
stress), the value of R is -1.0, which describes a fully reversed state 
of stress. The quantity Y , which is an explicit function of the crack 
length, is a geometric correction factor that accounts for the decrease in 
load bearing area during crack propagation. Ac such, Y can be satis- 
factorily appaoximated by a second order polynomial 
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y - aCa/t) 2 + e (IV. 7) 

where o,e ■ curve fit coefficients. 

Combining Equations (IV. 5) and IV. 7), and separating variables, yields 



Evaluation of the Model 

The fatigue life integral and its supporting equations were pro- 
grammed for solution on the CDC6600/Cyber 74 system. A trapezoidal inte- 
gration scheme was used to evaluate the definite integral in Equation 
(IV. 8). A listing of the computer program, FATLIF, is contained in Appendix E. 

Since prediction of fatigue life is currently accomplished by a 
stand-alone program, it was necessary to first exercise program "Bellow" 
to generate critical flow excitation velocities for a given bellows con- 
figuration. The essential input-output data for program Bellow is summar- 
ized in Table X. The reader will be able to identify the bellows and fl- -u 
input parameters that are common to the FATLIF program. Excitation velocities 
are shown for the first four modes. It should be noted that the fatigue life 
program accepts fluid pressure in psia rather than psig. 

To complete the input data for program FATLIF, it was necessary to 
specify numerical values for C, m, a, e, a^» and a c . The constants, a and 
e, were obtained by curve-fitting the correction factors for a single edge- 
notched strip that are presented in Table 4 of Reference 14. The results of 
the manipulation yielded 

a - 6.79 


e 


1.12 
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TABLE X 

INPUT/OUTPUT DATA FOR PROGRAM BELLOW 


BELLOWS PARAMETERS (INPUT) 


^lGMATCOfivOLUTE ‘ SIOThTTm) .158 

LAmBQA( Cuwv'OLUl E PITCH; JN) J.PbR 

~H(k£AN disc HEIGHT," IN) ~ ,3MO 

~T (C ON VUL DTE THICKNibSS/P LT, IN) ' ,00 8 

D’l ( InSI.OE *t)i AME TER, IN) H.ooci 

0 _ 0( OUT SIDE ~0I A.i£ TEW, IN) *.«28 

NCC'nUmUEW OF 'CONVOLUTES)" R.OOij 

NPLY ( NU^-itK OF PLIES) ‘*3.000 

eT young 1 s modulus; Lc/SQ. IN> P.fiHboE+O? 

~KA(OVEK'alL SPRING RATE", L6/ IN) 3 73 . 3 b b 

WHO h'Ch A T c R I'A L ' 0 1* NS I Tjr / C^CU .J_NJ> P.82 


f L UIO P AhAM E T Eft S (INPUT). 


pIpkessuh'e, psio'y 100.000 

Temp ( temperature; oeg n ?o..ooo_ 

" “KhOT ( F L 0*1 o JJc’N SITY," L3/CU.IN j ' 3 . b 1 11E -U i 

NF L U 1 0 ( 1 = G A s; 2 = L I Cl U I D ) . ** ~ _2 


Mode 

Hz 

THEORETICAL PERFORMANCE (OUTPUT) 

Flow Excitation Range, ft/sec 

Lower 

Critical 

u PP- e I 

i 

191 

7.0 

9.8 

16.2 

2 

365 

13.4 

18.7 

30.9 

3 

523 

19.2 

26.8 

44.3 

4 

666 

24.4 

34.1 

56.4 
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Specification of the values of C and o was impeded by the lack 
of basic crack propagation data in the open literature for the bellows 
materials of interest, i.e., Inco 718, 21-6-9 and 321. .As an alternative, 
for evaluation purposes only, the following values of C and m were ob- 
tained from Reference 15 for Type 316 stainless steel. 

C - 7 x 10~ 16 

m » 6.5 

The number of fatigue cycles needed to effect failure is strongtly 
affected by the magnitude of a^ and a c . In evaluating the model, the 
initial flaw size, a^, was chosen to be 0.001 inch. This value is be- 
lieved to be representative of a typical surface flaw. The crack length at 
failure, a c , was taken to be t/2, the validity limit of thp model. 

Based upon the above inputdat;., fatigue life predictions were made 
for the specific bellows geometry, fluid properties, and critical excitation 
velocities in Table X. The results are summarized in Table XI. For this 
example problem, the following observations can be made on the validity of 
the model. 

(1) The cyclic stress range, A o, increases with mode number 
because the product of flow excitation velocity and dynamic 
amplification factor is an increasing function of mode 
number . 

(2) The maximum bending strews is tensile at all mode numbers. 

The minimum stress is tensile initially but becomes com- 
pressive as mode number increases. In the presence of in- 
ternal pressure, a fully-reversed stress field is not 
achieved. 

(3) For this example, in which Type 316 stainless stee_ was 
employed, the maximum tensile stress in the first three 
modes did not exceed the material yield point of 42 ksi.(16) 

In the fourth mode, the maximum tensile stress exceeded the 
material yield point. 

(4) In this example, the model predicts high cycle fatigue when 
a^ and a c are 0.001 and 0.012 inch, respectively. 














At this point, realization of the full utility of this approach to 
fatigue modeling is impeded by: 

(1) The lack of basic crack propagation data for the three 
materials of interest. Currently, the fatigue data that 
are available from the materials manufacturers weru ob- 
tained using fully-reversed stress fields at room temper- 
ature. What is needed are crack propagation tests which 
yield crack growth rates as a function of stress intensity 
factor and mean stress over the range of temperatures of 
interest. 

(2) A correlation between the fatigue life as predicted by the 
crack propagation model adn experimental fatigue life of 
actual bellows in a common temperature and stress environ- 
ment . 

(3) The lack of an experimental definition of the flaw size, 
a^, that is needed to initiate and propaga' a crack and 
a c , the crack length at which failure occurs. 

VI . 2 Fatigue Curves 

Due to the limitations posed by the crack growth model, an alternate 
approach was developed to predict bellows life. Seven ordinate charts were 
developed (Figures 23, 24, and 25) from data listed in References 17 through 
27. The materials studied included 347 SS (a close substitution for 321 SS) , 
Alloy 21-6-9, and Inco 718. Data reviewed were mainly in the form of "cycles 
to failure" or S-N curves for various R values and for temperatures of 70 # F 
and -423 °F. 

Seven-Ordinate charts relate stress and stress ratios to cycle life. 
Most of the seven-ordinate data is based upon data banks maintained by the 
Department of Defense and the Federal Aviation Agency if its source is con- 
tained in the MIL-HDBK-5B. 

Seven-ordinate charts are convenient to use and they relate fatigue 
life in terms of mean stress which could be an important factor when predict- 
ing bellows life. Design or analysis parameters can be specified as stress 
amplitude, mean stress, maximum or minimum stress, cycle life, R-values, and 
A-values. (The A value is defined as the stress amplitude divided by the 
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3 



— — — = -423°F 

= 70°F 


Minimum Stress 
Maximum Stress 


FIGURE 23. SEVEN-ORDINATE CHART FOR ALLOY 21-6-9 
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INCONEL 718 
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FIGURE 24 


SEVEN-ORDINATE CHARI FOR INCONEL 718 
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347 STAINLESS STEEL 
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mean stress.) To determine the fatigue life, only three parameters are re- 
quired. These parameters are usually determined from the bellows stress 
analysis. 

The seven-ordinate fatigue data is built into the computer program 
listed in Appendix A. Each constant life cyci curve is modeled as a power 
law, or 

Cj “ B 0 alt B 

where Cj is a constant life value for j mean stress. Curves are devel- 
oped for mean stress levels of 0, 20, AO, 60, and 80 ksi. The alternating 
stress component, 0 a lt* * s i n the units of ksi. Simple linear interpola- 
tion may be used for intermediate values. 

The seven-ordinate curves are applicable for fatigue life predic- 
tions once the stress levels have been determined; however, stress indicator 
values may be used directly as a calculated alternating stress value with 
reasonable accuracy even though the stress indicator's intended use was to 
predict fatigue life with the aid of data presented in Figure 6. 
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BELLOWS FLOW-INDUCED VIBRATION COMPUTER PROGRAM 
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A.l Governing Equations 

The performance equations, which will be presented in this section, 
are based upon the derivations given in Reference 1. Therefore, detailed 
algebraic manipulations and derivations have been eliminated for clarity. 

Figure A-l illustrates a longitudinal cross-section of a typical bellows 
together with pertinent notation. The overall bellows spring rate is 

K A - D m E ^ (t/h) 3 (A-l) 

where E is the Young's modulus for bellows material and D m is the mean 
bellows diameter which is defined as 

D m - (D ± + D 0 )/2 (A-2) 

The elemental spring rate, K, is given by 

K « 2 N c K a (A- 3) 


The corresponding elemental metal mass of the bellows is 
“m * * Pm c N p D m [* a + (h - 2a)] 


(A-4) 


where p m is the metal density and the mean crown or convolute forming 
radius is 

a - (<J- t N p )/2 (A- 5) 

As the bellows vibrates in any one of its 2N C -1 longitudinal modes, fluid 
is accelerated within the convolutes. The process of moving the fluid is 
manifested as an apparent of added mass which must be taken into account in 
calculating the frequencies at which a fluid-elastic instability is likely 
to occur. This added mass is a function of the longitudinal mode number, N. 
That is 


/2 H c -1 -tT\ _ ( N-l \ 

°f “ m fl \ 2N c -2 j T “f 2 ^2 N c - 2) 

where 

ir p f D m h (2a-tN_) 

“fl * Ti 


* D m Pf h3 
m f2" 3 g 6 


(A-6) 


(A-7) 


(A-8) 


l 



l 



CROWN 


A-2 



N - NUMBER OF CONVOLUTIONS COUNTED 
C FROM THE OUTSIDE 

N - NUMBER OF PLYS 
P 

D m - MEAN BELLOWS DIAMETER 
m 

t * WALL THICKNESS (THICKNESS PER 
PLY IF MULTI-PLY) 

X • CONVOLUTE PITCH 

0 - CONVOLUTE WIDTH 

a • MEAN FORMING RADIUS 

h - MEAN DISC HEIGHT 


FIGURE A-l . BELLOWS NOMENCLATURE 
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In these expressions Pj is the fluid density, g is the gravitational 
acceleration constant, and 

<5 » a - 2t N p (A-9) 

The node number, N, ranges between 1 and 2N C -1. A reference fre- 
quency for a particular mode number can be defined as 

f o <« ■ (A - 10) 

The true modal frequency, fjj, is then obtained by multiplying the reference 
value by the dimensionless frequency corresponding to the desired mode 
number and system degree of freedom. Dimensionless frequencies can be calcu- 
lated as 


B i 



% 

; i* 1,2,3,.. .2N C -1 (A-ll) 


Alternately, for purposes of hand calculations, the dimensionless frequency 
factors may be determined from Table A-I. 

It has beer* observed that flow excitation of a particular mode can 
occur over a broad range of fluid velocities, which is termed the "lock- in- 
range." In fact, if the modal frequencies are sufficiently close together, 
the lock-in ranges may overlap, thus producing nearly continuous excitation 
of the bellows. These lock-in ranges are estimated as follows. Extensive 
experimental studies have revealed that the Strouhal number provides an ex- 
cellent means of correlating the vibration frequency, fluid velocity and 
bellows geometry as shown in Figure A-2. The Strouhal number is based on 
convolute pitch, O. For a bellows having a convolute pitch-to-convolute 
tip width ratio of X/a, three values of the Strouhal number are indicated. 
Peak bellows excitation corresponds to the curve marked B a cr i t f rom which 
the critical flow velocity may be calculated, i.e 


v 00 

crit 


f N g 

S 0 crit 


(A-12) 


Similarly, the upper and lower values of velocity, which define the lock- 
in-range are obtained from 


v (N) , M 

upper S C£ 

£t£ 


(A-13) 


and 


v 00 

lower 


(A-14) 
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FIGURE A-2. COMPOSITE OF ALL STROUHAL NUMBER CORRELATION DATA 
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The stress indicator is a relative measure of the stress intensity. 
Two methods of calculation are allowed in the computer program. The first 
method, and the more exacting one. Involves a greater number of calculations 
and a substantial amount of input data. The second method incorporates a 
"Universal CpQ Function" and, due to its data compression requirement, it is 
by nature a more conservative calculation, i.e., the SI values will be high. 
These calculation methods are given as: 

Method I: Conventional Stress Indicator 


SI - 




C f C e P d 


N r 


- (h/tV 


(A-15) 


where * vortex force coefficient which is a function of \/o and is 

obtained from Figure A-3. 

C e * elbow factor to account for above average forces exerted on 
bellows convolutes if an elbow located immediately upstream 
of the bellows. 


* fluid dynamic pressure. 

Q * dynamic amplification factor. 


The bracketed term in Equation (A-15) is termed the "bellows opera- 
tional parameter". This parameter is used in conjunction with the bellows 
specific spring rate and Table A-II to determine the dynamic amplification 
factor (Figure A-4), where specific spring rate is defined as 


k a n c 


Dm N p 


(A-16) 


The computer program currently calculates the stress indicator corresponding 
to the critical flow velo ity defined by Equation (A-12). 

If the internal medium is a gas, a radial acoustic resonance condition 
is likely to occur, wherein the acoustic pressure fluctuations couple with 
the vortex shedding process to produce a force amplification that is signifi- 
cantly larger than .jould be predicted by the value of Q obtained from 
Figure A-4. Physically, these pressure fluctuations are attenuated at approxi- 
mately a constant rate for all vortex shedding frequencies less than the 
radial acoustic resonance or cutoff frequency. In the vicinity of the cutoff 
frequency, the increased amplification must be taken into account since it 
results in much higher bellows stress levels. To this end, the first mode 
radial acoustic resonant frequency is obtained from Figure 5 for a particular 
bellows geometry. This cutoff frequency is then compared with the predicted 
longitudinal modal frequencies. The predicted Q value from Figure A-4 is 
modified by a suitable constant for all longitudinal frequencies that exceed 
the cutoff frequency. In other words, this adjustment of Q states that the 
radial acoustic resonance is capable of coupling with higher longitudinal 
modes not just at the condition where the frequencies coincide. Figure A- 5 is 
valid for convolute pitch-to-tip width ratios of 1.4 to 2.0. These values 




FIGURE A-3. SUMMARY OF BELLOWS VORTEX FORCE COEFFICIENT EXPERIMENDAL DATA 
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TABLE A-II. APPLICATIONS INFORMATION FOR USE WITH 
Q VALUES DATA IN FIGURE A-4 


Specific Spring 

Number 

Internal Media 

Cuive 

Rate 

Plies 

(see Note 1) 

No. 

All Ranges 

1 

low pressure gases 

1 

over 2000 lb/in^ 

1 

high pressure gases, light liquids 

1 

over 2000 

1 

water, dense liquids 

2 

under 2000 

1 

high pressure gases, light liquids 

2 

under 2000 

1 

water, dense liquids 

3 

over 3000 

2 

All 

3 

2000 - 3000 

2 

all pressure gases 

4 

under 2000 

2 

all pressure gases 

5 

2000 - 3000 

2 

all liquids 

5 

under 2000 

2 

! 

all liquids 

6 

over 3000 

3 

All 

4 

2000 - 3000 

3 

All 

5 

under 2000 

3 

all pressure gases 

5 

under 2000 

3 

all liquids 

6 


NOTE 1: Low pressure gases will be defined here ac being 

those gases below 150 psia. Light liquids will 
be defined as having a specific gravity of less 
than 0.2. 



DYNAMIC AMPLIFICATION FACTOR 
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FIGURE A-4. DYNAMIC AMPLIFICATION FACTORS FOR 
VARIOUS BELLOWS APPLICATIONS 
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correspond to total convolute thickness of 0.3a and 0.0a (theoretical zero 
wall thickness). In addition. Figure A- 5 is valid for fluid danping numbers, 
1)^, of the order of 10“^ where D n - v / r i c o» v - fluid kinematic viscosity 
and c Q - isentropic speed of sound. 

Method II: Calculation of SI with Cw* Function 


Calculation of the stress indicator may be greatly streamlined if the 


universal Cp function shown in Figure A-6 is incorporated as follows: 


SI - C F * 


N c \ 
n n ; 


C e (h/t) 2 P d 


(A-17) 


Note that the calculation requires the use of only one curve, and hence, 
this method is favored for hand calculations; however, if JCFQ is set to 0, 
the calculation is performs: by the computer code. Input cards 9 through 15 
may be blank cards. 

Calculation of fatigue life is accomplished in a subroutine called 
XLIFE where the input parameters of material type, alternating stress, and 
mean stress are manipulated in conjunction with a "Seven-Ordinate" fatigue 
cha^t to determine the bellows expected life. 

The current version of the program assumes a mean stress of 0 psi; 
however, several simple program statements could be included to account for 
internal pressure and slight angulation. Room temperature conditions are 
assumed, but these conditions predict shorter life expectancies than cryogenic 
conditions. 

The room temperature conditions compensate somewhat for unknown work 
hardening effects. From the limited amount of data available (AFRPL-TR-68-22) , 
it is generally shown that hyd~oformed bellows life expectancy is shorter 
by one order of magnitude than that of a coupon made of the same material. 
Therefore, it is not advisable to expect longer bellows life due to low 
temperature operation. 

A tvpiral Seven Ordinate Chart is shown in Figure A-7. The alternating 
and mean stress ordinates are used exclusively in the bellows code. Each 
constant life curve is represented by a simple power law of the form 


Cycles * B a a ^ m (A-18) 

B and m values are obt Ined from data cards 16 through 21. For example, 
the cycles to failure for INCONEL 718 operating with a mean stress of 0 psi 
at room temperature is 

- -5.1097 

Cycles ■ 2.1410 x 10* a a lt 


Similar curves are generated for mean stress levels of 20, 40, 60, and 80 KSI. 
A linear interpolation process is used to compute cycle values between succes- 
sive 20 KSI mean stress levels. 
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FIGURE A-7. SEVEN-ORDIN/TE CHART FOR INCONEL 718 
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A. 2 Equivalence of Theoretical and Computer Program Variables 

This section is intended to establish the correspondence between the 
analysis variables presented in the previous section and the computer coded 
variables. Internally geneiated variables as well as curve fit coefficients 
will be discussed in subsequent sections. 


Analysis 

Nc 

N P 

a 

X 

h 

t 

Di 

C o 

E 


N 

X* 

a 

C e 

Pf 

Dm 

V 

a 

®m 

mfi 

m f2 


m f 

<5 


Computer 

NC 

NPLY 

SIGMA 

LAMBDA 

H 

T 

DI 

DO 

E 

RHOM 

KA 

CE 

RHOF 

DMEAN 

K 

A 

MMEYA7 

MFLUID1 

MFLUID2 

MFLUID 

DELTA 


Comment 

Number of bellows convolutes 
Number of ply? 

Convolute width 

Distance between adjacent con- 
volute crowns 

Mean convolute disc height 

Thickness per convolute ply 

Bellows inside diameter 

Bellows outside diameter 

Young's modulus of bellows 
material 

’’allows material density 

Overall bellows spring rate 

Elbow loss factor 

Fluid density 

Mean bellows diameter 

Elemental spring rate 

Mean convolute forming radius 

Elemental metal mass 

Apparent fluid mass at low mode 
numbers 

Aoparent fluid mass at higher 
mode numbers 

Apparent fluid mass 

Internal convolute width 





Analysis Computer Comment 


S °l 

STLO 

Strouhal number defining the 
lower and upper bounds on 

S °u 

STOP 

lock-in-range 

S °crit 

STCRIT 

Strouhal number for severe 
excitation 

V (N) 
lower 

VCMODE.l) 

Lower velocity bound on 
lock-in-range 

v <N) 
crit 

V(MODE, 2) 

Flow velocity for maximum 
excitation 

V (N) 
upper 

V(MODE, 3) 

Upper velocity bound on 
lock- in-range 

Cf 

CF 

Vortex force coefficient 

C F * 

CFSTAR 

Envelope stress coefficient 

SSR 

SSR 

Specific spring rate 

Q 

Q 

Dynamic amplification factor 

SI 

SI 

Stress indicator 

°alt 

ALTSTR 

Alternating stress 

a m 

MEANSTR 

Mean stress 

w co r i /c o 

FNCO 

Frequency number for first 
mode radial acoustic resonance 

w co 

FREQCO 

Angular cutoff frequency for 
first mode radial acoustic 



resonance 

Co 

CO 

Isentropic speed of sound 

Y 

GAMMA 

Ratio of gas specit ic heats 


A. 3 Curve Fit Requirements 

When predicting the performance of complex systems, it is frequently 
necessary to describe experimentally observed relationships between two or 
more variables through the use of empirical expressions, i.e., curve fits. 
In predicting bellows flow- induced vibrations, it was necessary to curve 
fit the data shown in Figures A-2, A-3, and A-4. To this end, all data in 
these figures were fitted to a hyperbolic equation of the form 
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y » — ~ — + b + dx (A-19) 

X ** 3 

where k, a, b, and d are the coefficients to be determined. Coordinate 
pairs are input to the fitting routine, and the resulting equations are 
solved simultaneously for the unknown coefficients. A listing of the curve 
fit routine is included in the next section. Note that there is an option 
for either a four- or eight-point fit. It was necessary to use an eight- 
point fit only for the cuives labeled 1, 2, and 3 in Figure A-4 (Q-surface) . 
Curve fit coefficients are supplied on input cards 6 through 15. 

A. 4 Compute r Program Structure and Listing 

The computer programs listed in this section were written in FORTRAN IV 
language. Ir. the form presented here, the programs must be compiled each 
time they are submitted to the computer; however, multiple runs can be 
accomplished at each submittal. The user of this program may find it more 
convenient to compile and store the program on tape, thus necessitating minor 
program modifications. 

Four program listings are contained in this section: 

(1) MAIN (PROGRAM BELLOW) 

(2) Curve generating routine (CURVE) 

(3) First mode acoustic response frequency (ACOURES) 

(4) Fatigue life routine (XLIFE) 

The source deck for the performance program consists of a main program in 
which a majority of the calculations are performed and three subroutines: 
CURVE, which is called from the main program, and it contains the logic for 
selecting the appropriate curve on the Q-surface (Figure A-4); ACOURES, which 
evaluates the first mode cutoff or acoustic resonance frequency as a func- 
tion of bellows geometry, and XLIFE, which calculates the bellows life ex- 
pectancy based upon seven ordinate fatigue data. 

The execution structure of the program consists of the following items 
in the order presented. 


(1) 

Program control cards - number and type of 
varies with the user facility. 

these cards 

(2) 

Main program designated Program Bellow. 



(3) 

Subroutine CURVE 



(4) 

Subroutine ACOURES 



(5) 

Subroutine XLIFE 



(6) 

End of record (EOR) card; multi-punch 7-8-9 

in 

column 

(7) 

Data package containing one or more runs. 



(3) 

End of file (EOF) card; multi-punch 6-7-8-9 

in 

column 




PRuGRAM 8ELL0*( input, OUTPUT, TaPEsOsInPuT) 

THIS PROGRAM GENERATES * THEORETICAL PREDICTION OF THE NATURAL 
VELOCITIES NHlCn PRODUCE FuO— INOUCED VIBRaTIONS(EXCITaTIOn) OF 

-PREOUENCIES' FOR'”a^I V En — 6'E CUD hS — In CUID"InG — TRE — FLUID - FUD* 

Tn£ _ 6 £ LCD *3 '~N ~ATURTrTQNGlTgDTNAL~NggnT ' 


IRPUT 


■^FLT<r«nrc^cgLTT6 - ^r^n7srTxP6TnNTWiTiTi^TtTR>ffW'irA)Tiu — 
i s r fiE~ffVE r rare ecto ^“swing 'Ratet ’ ls'/In 

— NFLuTO - *~C( nxsm (CTSDlffT 

~ Nti EG' =-KU M6£ R“OF“B E CTO^-iTjn (TnWr^LTrDRrrsnjF^EI (fO^^NCri — 

J hTT^ _ n 0 hb£R QF' CuRVES RCcESSaRt to DESCRIBE 0 SuRFaCE 

NC * nuhPER' OF BEll BnS CQNvQLU TES ' 

~nPU?"«~nun¥E ^~ 0F~>lVS in ThT~ BELLOWS cqnvqlutes ■ 

^sigha » conVoTutE »IQTh , In. 

LAlig D A_ «_ DI ST ancE~B&ThTEN aPJaCEnT CONVbLUfg CRQXN3> Tn. 

_h _* _n£ a N_OI sc H EtGftT, in. * 

t * thickness p'e'r convolut e ply, in, 

Dl * BELLO »'S'TnS7oE - 5 I a HE'T l R , In. 

J>0_« e'ELLO» S~~QUTSldE DIaheteR# IN. 

E~* fOUNG's HODuLUS OF T HE SS L LOwS M ATER IA L, LB/SO IN. 

:*»ho>t« hEIGh T DENSITY QF'ThE SELL OKS haTV'RIAL, LB~/CU IN. 

c£’ * disensionleSSTElbo- FALfOR 

IF* nf'LLIO * 1 (GAS Th£ FWTcT GAfEOUAflON OF SfTfTTS USED FOR' 
CaCColaTIng GaS OTnSiTtftT Th-. stTte OIfinEo by p 'and Temp*. 

" IT 1 S VSSu HE'D "thaT^ThE G aS — VRqp I (f TTe S are knJThn a T“a~ REFERENCE 

S TaTE D fF I ' n£F ~B r rhOEr1f7~>'REF, and TR~S~f . 

P r- i — G a S ~PH E~S~S~UR£, ^S lG 

TE’mR z GaTHYhFeiTaTu&£, OEG. F. 

— pref-awtref— ^FER r^r^*s'' , 'stATE;“Ps'!A'ANbn aTtrrr. 

RhuF R£ FTTS a S" DENSITY a T RTF' E'fllfN c TATE, Lb / C U FT. 

— gAMHA---r^iTiD^F^prcrrFTr^rA'Ts-ROfnrA's 

“tf-mflui o "* ■ s cn wuioyr t he crcgrr'bT.yr! ty-wst-pf'knE'-n ~ a rwto r rrr 
The' LI uu 1 6 tfTJTfO* - anD T £mF5 

— t>— ■ “L raonr^RTssuRr. rtstg 

T£*P * l roul"D T I'nPTS a T U R E , OEG. F . 

R nOF-r TI 9L< TD - DETJ5XTT A'T~ F - a"ND — TE m FTTE7CITF7 7 

HTLTia. T E R I A iTTN 0 TCATO RTITT nCC T^TS ,~ £TtCCO V~8 r^T-q~3T3gVSS7 

- TTU B K7‘ST'URAVSTiJPBT5TgRD'~* — CURVE - FIT - COEFFICIENTS' PgR - jPPER''adUND' 

ON" 5TR0UHAC - WONBER - FS7""L* tf 3T>A7STSRA 

-^LT.R7STLXA-,^TC06T5TCap-^~SPT e ~ rs ~ ie DVE-EYCE' PT ~ L . u VEF~B~0tTOD 

-^ T CSTTK;STCRTriT5TCFTTB75TCRTT5 - «-3'ARE - TS - A'BT)VE - FXCEPT-F&R - 0PTlHyM- 

O R LRr ncmTROOTAX -WPB'prTOW" BECU0F5 - EVCTTITItfN 

~TFiC7CF~A7CF37i:rD— fTORVE~T i T-XgE > F7nTUTr'FOTr-VORTErrORCr 

COEFFICIENT 


-OJa 0 ) ; 5 a CT)-, OBT JTTBD O rv~CO RVE" 

AHPLTFICATTDsr FTCTOR CO' “SURF ACE 

->n«— DI NENSIONLES S — i aT 0 BTC - FR WUEnCVaS a FONCTniN.- 
OF - **OD£- nuhBeR'FoT nOIG^TClOhTTOngITUU'i'nIL 
DEGREES’ Or’FPEEDOKT; 


i’HE — CTTNAMIC” 


>Fva;NTLi-«--rRirTn'^S70facr^rTRTrTCg»nTiwi>icr-vTO3ES-or-N tn 

— C AL'CUCaTINGTTFE “CYCCE S 7~P TC~I 3~> r rTE RT1C - T N 01 CTT0R7 

■> (tt-tl -) - f TworoTHERrrcNfL hatrit-coftaTn nre-v atuettifti n 

— OTX UL'aTTMGTT rFr~CTC LlS T~ > T L ~ I S~VaTI RTaC - IWaXaT O'RT 

■S1TBR0VTT NE~» CTF ET^VCOCTTES - rHE^JMSTER-OF-PREDTC TTlTCI FE~ CYC LE S 

— G rVE ^ - At TCRNAnNG~'ST»?rSSrvS’I )~7~REATrSTROSTKTTr, Th7an5 ~B~ Talu£S 


OMQMAL PACE K 
Of POOR QUAL.TY 
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000003 

0C0003 

00U003 

DIMENSION FAEQ(35), V(3S,3),3I(3S),TFAIU3S) 
DIMENSION Qn(b).QACb),66(b),QD(b) 

DIMENSION XM(5.3),6(S,3) 

oWool 

DImEn'sIOn TlTuErO 

oooool 

DIMENSION xNHfsjTJftnTsl —— 

oTJb'5o3 

»F*L n£$lo«(D » mEsTTI (T) , m£aNI.(SE*T “ 

oOoTTcfl 

REAi. NC#NFCV»L*MSDT7SFirUTirrrMFtUl64#M^l.UI0»M«ET4l.»l<A,K,M*iS 

OdtlOdS 

Tf£AL *EANST»l 

000003 

CommCn n^TTS^R, NP fcTlDTdTR'Hd*' 1 * JcurVE, &5P.Q, CTCLE , MRET 

030003 

&ru MEsLdN/br»sELO',*Mb :oT*Mt*oT7 

OUuOfil 

OaTa MESSr/*MAdO'»#*r!g io»*h6*oT7 — - 

oooool 

iTaTa nEAr4LC74HHE*KT,4H $Tfi#bH£$s ,4HHl 7 ^ “““ “ 

c * • » » » * • •**•«****•**•«••* 

000003 

1 R£aO 1000»TITlE 

oinioTT 

fF(£0? , bO)S, 10 

00001* 

S STUP 

00001 b 

10 REA'D 1030, jFun5,NFLOtd,NDE5, k 'N*X, JCFOfMfc “ 

00003b 

READ liai)#Nd7N>L?,3l6MA,tAM80A,H,t 

00005b 

READ 1010,01. DO, E.RmOm.ka.CE 

00007b 

60 TO (11,13), nFlUID 

000104 

1) READ iOid,F,TE«R,RREE,TRE^,RMO^RE^,CA!**iA 

000124 

60 TO 13 

000135 

13 RE*D 1 Oil), P, TEwP, PrtdF 

00013? 

13 »E*0 iU80.3TU»K\STU(^7sfoPFrsfyP0 

0001ST~ 

READ 1 0 80 , STUOkTSTlQA, ST lda, St LOD 

oBffTbl 

READ i08d,S>CRrTA,STePITA,STePfTd,STdftlTD ~ 

000203 

RE*0 1080, CFK,CFA,CFd, CFO 

0l031> 

READ 1080, (0<cT),6AiJ),ffo(J),00(Ji, J*1,J*AX) 

000240 

IS READ 3dbl,(XM(I,I),)«l#sT " 

000252 

READ 3003, (6(1,1), 1*1,5) 

o3o7b» 

REaD 3001, (XH (1,3), 1*1, 5) 

0002?b 

REaD 30037(1(1,8). 1*1,5) ” 

0003X0 

3*£ac 30di,(XM(i,3),l*i,s!) “ ' 

000322 

READ 3cd3,(DTt,3).I*i.S) 

c * * Calculation of natural frequencies and excitation velocities * • * 

00033* 

25 PI*3.1415R27 

00033b 

G*32. 174049 

00033? 

0m£an=(01*00)/3. 

0003*3 

60 fO (30,35), JFtTS ' 

0OCTS3 

3'C AA*DhrAN«E«lNPCT/NC)*(T/H)»»S 

ootnSs 

~ TS k«2*»nC*ka*i2. 

0003b* 

a«iSIGma-T*aPlY)/2* 

OfiOSbb 

^StTAL#Pl*RH0NTt«NPLlr»DN6AN* r &I*A*H-2.*A'i /& 

odds?? 

60 TO 1 3b , 3?) , NFCUlO 

000*05 

3b R*tOF*RHOFBEF*(P*l*.1>)*((TREF**bO.)/(TENP**bO.))/(PPEF«i738.) 

000*1 b 

GO TO as 

000*1? 

37 RnUFxNHQF/ 1738 . 

000421 

38 HFLUl5I*el*«H0F*6NEAN*H*(2.*A-T*NPLr5/C2**G) 

000*33 

DE£^a~*SIGha-2.*Y *nply 

Odd* 3b 

MFuOl D3»PI*Ri*dF«0*N£AN* (m»»3) / ( 5 . »6»0£lTA) 

000 * »s 

x*LA<46Ca/s!6ma 

00U4V> 

5 t i. C « s t L 0 K / ( X .1 f lo AT^5TO5f7XTTtf5^ 

oTo*T? 

^Typ*sTuPK/tx-sT0PAT6ITuA57rTyPD»x 

5TT5T51 

STD R i T M s T C Rl T k 7 (x - 5 T CTTTaT + 3 TcRTTB~*TfcFn b • x 

000471 

ahOoe*I.O 

000*?S 

DO bO MOOE»l,NOTtf ’ 

0004^4 

0E«rl«3.»NC-3. 





000*77 

000500 

ooosei 

000512 

IF(0E*Fun£.0) 60 TO 3* 

OE«Ft.si. 

3* "tFUOlD*(MFLUl0l«(2.*NC-l.-»H0DE)*i*fLUl02«CAM03E-i.))/0EHFL 
*A5 SshFi.UI0***ETal 

fnren* 

8 1 « $6ft T ( a . * Cl . tfCsrOl* ( 2 . *Nt-M(JOE ) T7TT~. **t ) ) ) ) 

owm 

F RE8X MOOE ) *'S6a TOO SX5IT*&T7T3 . *T f T3 

ooms*3 

DO 35*^81,3 

OioS*3 

eo t6<*o,*s7IoT7J “ 

BW5T5T 

*3 vc«OOE» J)*FaE6( 1*01)13 »fI5* Ti ( i 71 

0005*0 ‘ 

' 60 T O’ 55 

0005b 1 “ 

*5 9 1 i*Od£ . Jf) *F»li3XR05S^»SifeHA/CStCRIT«l2.) 

o0o57l “ 

to ro 'ss 

" O0o57i 

50 V ( nOOE ,J ) xF»£S( h'CoE) *3 IG* A/ ( STL0*12.) ' “ ' 

oootfli 

55”C0nT InUE 

tiQUblTJ 

— sir ‘AiGDEsAftooevn 

£»»•«•««****•»«•,•,»***•*•*•«*•*•**« 

C T'm£0k£T2CaI fffHESS^lNOTJrrtO^TT^P fTcTL* ST&Q uhal N u*0£ff 

OOl'b 1 U 

CF«CFK/(X.CF>J^CFB^:Fa>X 

OCOblb 

ssasR*»nc/(6f*EAN*N?T , n ! " 

” OtJnbJI 

C*lL curve 

OOtib? 2 

IF'CnFluID.ES.I) *5 , 70 

OQOb27 

b5 AI8DI/2. 

OQUbSl 

HtflSh/SI 

00 U b 3 3 

C0*56R T (6A'*MA«(R»l*.7)«6/(RM6F*ia.)) 

000b** 

C*LL *C0URES(nSIfRI>CO»FREOC>'',O*CJUS T ) 

000b* 7 

'JO A«00£*i . 0 


OOObSi 

50 T ( ?3 # 1 bS ) p JCFQ 

OOQbS? 

>1 oo iso *ode*i»nTeg 

OODbbi 

PA Ft A* f*Q 

000b b* 

IF(PARA.Gf..O?S> GO TO *0 

OOlibb 7 

CPSnSsi.o 

000b 70 

30 TO 1*0 

000b7r> 

*0 IF(PARA.GT..lb) 60 to 100 

— OlJTJ¥T* 

CFST.*.»x.b*?S*3* 1, *R»'»C-.13A0Sb ) 

oooloo 

GOtS 1*0 

ffffoToo 

ico iF(p*fit*.6f773) g 3 to iio ■ 

ouoTo^ 

C ? S T 4 n 8 , 2 7 5 cj i b • ?"aR * »TriTb*S 0 2 ) 

oTTPTiTi 

Gtf TC r*o 

tFBF7T5 

ITT; 1 rc p « * . c T . . 7j GO TCT liSo 

60071^ 

CFST4?*7.3S33* , f*P'iRTi*J.08l*S? 

55TT7TD 

60~T0 1*0 

” mrm * 

FTP AR A; GT72T3 G"D"'T 0 Tfo 

UTOTTT?^ 

wFStTPs^7^0lib7*PT^A*V7nT?T> 

' 3WTJ5 

GO~T?r 1*0 

oBZFTyti ” 

lie CFST *R*b # U 

oxrcni? 

ffo Ba«8CE.Rn0FT(V(M00ET2) ** 8'f* (1 ( H7TT* «Y) 7Tf77T e .~*PL> • Gf ) * fCF 3 T * ») 

WHtt 

_ 51C«00E) xft'ON* (NCTa'^OC'E') 

0007S1 

• CFii*cFST aR»nC/ Afi5i)E 

* otr&T?* 

PRINT 1*00, CFO 

OfTQTbl 

1*U0 FoffKATtSH CF0*,Ci2.b) ‘ ~ 

0007bl 

ahOOE*AHOOC*1. 

oBoTbl 

ITo Continue 

ffaUTbb 

63 TO 205 

5MTbb~ 

TFs 00 2o 6 KOOfmTTNTjEC “ ~ ' 

cfftTmi 

eOP*CF«CE*R'<OF«iV(tiCnrE, 4)**4)«((H7T)»*i)*i2./(i. »npi.t*6) 

odioo* 

OxG* ( JCUflVO / fSOP^ff* O'fcuR VED 7SB(JtjfiVe3*S0A*5&tTarRVE 3 

GoiiilS 

IFT*Fl'uidTEQ.Y)T?S» r#is 
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00X017 

00102* 

00102 ? 

0010 31 

dOi 633 
6om»T~ 

~ " ofliiuf 
66 i 111 

001115 " 

OOl llb 

offllfo 


001122 


00112 * 

00112V 


00113b 

~6STl*b 

00117 0 
0011 7'l 

*60121 3" 
00121 * 



0012bl 


0012b* 
00 1270 

opisoo 


1?5 If (PREGCHODEVGE.PPEOCO) Q«Q«OiOJUST 
180 SUhOPE)«BOP*C 
*POOC**NOOE*1. 

20C CONTINUE 

"3nS" PRI N T~2 o'uoT^n T Cl ~ 

PftTKT • 1 d* j;s I C^uiMloTTrCTTO' ITMTHTTSPI y7? 7* A*; * hO*7p7t e hPT* HOPVftr 
TCuIO.pTl 

PPI N t* 10 SO — — 

MJTn f 1 6 b o ' ■ - ~ ~ 

nEJnsTP«0V5 

00- PV H00E*1*ND£G ' — 

*CT ST ft • s i (hoPI i / iooS C 


00 101 1*1,5 


xan(n«XN(I, PTL) 

e»N c i >*sri7NfLi 

TnT continue 

C*ll xuT f ET5tft«, S phTACT 3~ttt , heanstr ) 

GO" To CIVIC, 1 300 , lJtOo, 210T7 n0TT 

„ flo p»Ist" i 0 ?o7m6'6e75TTh06e ) « f fteo C * ooe) , v ( hooe , i ) , v ( hope , 2 ) , v c nooeTTjT 
$ - c'TCtE 

GO fO Pft — - 

_» 3HQ~ ~E"IV t 13 ns, M OD E, SI (NOPE) ,pp fe6(MQ0E) , v ( hO 0E , l ) , v ( hooE , 2 ) , wfi wOOETTST 

S pfS'LOi* ‘ , 

GC TO ** " ‘ ~ 

l^OC Pf»lNT~l305,H60E~,SI(wOPgT7TQro(HQOl),V(HODl,l), v c H ool 72777c mod £77)7 ■ 


PESHI 


GO TO ** 

1*10 P*lM i * 1 5 , H ODETS M H CO E T , F 0 E 3(*0DE) , V<H(jOE, ; 


PE4NLS 


CONTINUE ' 


IfLCNfwUjO.Eo. 1)215,1 
215 P RINT lU*U,fftEQC0.04PJUST 




P0NP4T(8*10) 


001301 1010 FOSP*T ( b£l2. b) 

0013C1 1 020 FORMA T ( IOFV.Xj 

0013 01 1 03 U_F Q.SH* T ( b 1 3 . 

001301 1 0*0 fO» "XT ( 1 NO ,2 * X , 1 8 - 9El.L0»S P*R*m£T 

% ”~_'lHO • 1 * K,2P~»*Sl6N*( CONVBtUTg blQTH, IN) , 1 1 X , Fb . 3 , / 


1 2 ?nC*N 8 Qt (CON VOLUTE PITCH, lN),inx,Fb .3,/ 

1 P »,23'hh(hExn PIS . HE IGHT, IN),l *X,?b 

IPX. JOhTC CONVOLUTE thT£Hn1ss/pi.t 


N),T3»,f».3,/ 












A-21 


S 3b*,2HnZ,lHXf SHLOwES^XrlHCRITICAL^X.SMUPPEXrbX^MSCC) 

001301 lObCJ FO«*»T(lrlO) 

001301 lO’O F0*w*TC3X,I3,8X,Ell.b,bX,Fq.3,5x,3F11.3>bX>Eli.'O 

001301 1080 FO«* ! ATC*El't.S) 

~o~u~ n ^I' — M Fr~FQ*"i r( 5f 1-P7 35 

“WI T H Io4(T \ * * » T (7 /77 3 * ,TbWT 0 IT IxiTTfO JSTTTTfn 

~~~ i . a # § *0 * o'j Cist » , F t. fj ~ 

bOltui n u5~F o«hVt 7 YxTr?7Tx7Eil.*.^x.Fb.t,sx, 3Fii.3>bx,3**) 
o6i3oi r*OTFO»H* TflxTl 27*771 ii.*, <*72 OTtx.ff i'l.S.bx, **-o 

001301 Th(5 


41 N~f~ FftrOUlNfT^ 
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000002 

000002 

SuftROuTlNE curve 

01 PENSION DIHFRE0ftS>.FRE0(25),V{3l, J),3IC»5) 
0 I"En3I0n 0K(b),uA(b),Q8(b),UD(b) 


flit) 

COMMON'* RCTTS3 RTSFlJT 07*7* hDTTJCUR^ETB'O R7TT,"CY£LE7*tfFr 


OOGOD2 

pe*c nc ; npc YTir-iffDATWF LuroTTHR cinoTr^iTnoTWHrr a l. kitk 


oTooii? 

IF (nply. EOTTTTI io f io 


mnnnr? 

TiTTF (KPL Y". E 0 .■J7yr(572B 


— 5BMTf — 

SCI F C 3 3 R . G T . 3tfB U?y 3TPF0 


oBooT! 

SoTCuRvE** 


■mm 

atTUKN 


0TTO1 

*5 IFT2O00. .LE.ySR. AND.35R.LE.3tFoO.)5fl*bO 


00003b 

50 JCuRVE«S 


0U0O3? 

RETURN 


0000*0 

bO lFlSSfi.C?. 2000. .ANO.NFLUlO.E®. 1)70*80 — 


oOoosi 

” >0 JCURvE«5 


000 ns? 

return 


offooTS 

80 JCufc Vt*b 


oooos* 

Return 

. 

ooooss 

"C IF(SSR.GT.3000.)100.110 


000003 

IOC JCUR vE*3 


0000b* 

return 


OOOObS 

no 1f(nfl'ui6.eo".i)i4o,iSo ~ " 


» HaOB 



000 10*> 

130 JCUR vE*4 


0U010S 

RETURN 


0001 0 b 

1*0 JCuRvE*5 


00010? 

RETURN 


ooouo 

lSu I F ft 00 * 0 . . LETfslTTNinTJirTr. 300 0 . ) l bO , 1 70 


oooiaa 

IbO JCu«vE«S 


000123 

RETURN 


000l?4 

1?0 JCURvf *b 


oooies 

RETURN 


00012b 

180 iFtNFLUlO.EO.l.ANb.P.LT . ISO . ) 1 *0 , 2UQ ~ ” 


000140 

ISO JCURVEsl 


0001*1 

RETURN 


0001*2 

200 fF(SSR.GT.2000.)2l0,2b0 


000150 

210 00 TO C220ft30),NFLUID 


OOOlSb 



offcTIT? 



OOOlbO 

& 

9 


OOl'lbS 



b JOibT 

2*0 JCUSvE«l 


otriTiTo 

RE t uRN 


oWTTI 

25o~JCuRvE*2 


























X 



R 

-> 






SUBROUTINE *COURES(X, Y,Z,FRE0C0,9*0JUST) 

? 


000010 


, 

; 

oooon 

MR I *X 

1 


' OObdl* 








1 
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00000? 

0011007 

—ffffwr 


SUBROUTINE XiIFE(XN'B'ALTSTR'MEANSTR) 

0INENSI0N xn(S) >B(S) 

CU"*0N nRlV»SSR.nflUID«R»RnOF. JCU»VE»B0R#Q»CTCLE»HRET 
"REaTTnEanSTR ' 



OOOObO 


OOOObl 


0000b? 


OOOU 70 


000071 


000072 


000073 


00007b 


$.10 HR£T»b 


RETURN 


£NC 
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A. 5 Data Input Package 


Instructions for preparation of a data Input package are located at the 
beginning of the PROGRAM BELLOW listing. An experienced programmer will have 
no difficulty In constructing the Input, but for the Inexperienced user the 
following supplementary remarks may be useful. 

Input Card 1 


This card is an identification card on which the user can place 
information that will aid in Identifying and classifying the run. 
Any alpha-numeric characters can be placed in columns 2 through 
80. Column 1 must either contain a 1 for printer carriage control 
or be left blank. 


Input Card 2 


Word 1 
Word 2 
Word 3 
Word 4 


Word 5 
Word 6 


JFLAG 

NFLU1D 

NDEG 

JMAX 


JCFQ - 1 

- 2 


MTL 


See Program Listing 
See Program Listing 
See Program Listing 

Is the number of individual curves necessary 
to describe the Q-surface (Figure A-4) . As shown 
in that figure, JMAX - 6. If future data indi- 
cate that more than six curves are necessary, 
then the dimension statement pertaining to Q 
must be altered accordingly. 

(Use Method I Stress Indicator Calculation) 

(Use Method II Stress Indicator Calculation) 

See Program Listing 


Input Card 3 

Word 1 NC 

Word 2 NPLY 

Word 3 SIGMA 

Word 4 LAMBDA 

Word 5 H 

Word 6 T 


- See Program Listing 

* M 

^ II 

« H 

It 


Input Card 4 


Word 

1 

DI 

- o a Program 

Listing 

Word 

2 

DO 

ft 


Word 

3 

E 

_ tt 


Word 

4 

RHOM 

_ If 


Word 

5 

KA 

- May be left 

blank if Jr 

Word 

6 

CE - 

1.0 - See Program 

I is ting 


Input Card 3 

Wov^ 1 D 

Word 2 TEMP 

Word 3 PREF or 

RHOF 

Word 4 TREF 

Word 5 RHOFREF 

Word 6 GAMMA 


See Program Listing 

t» 

It 

II 

II 

ft 


- 1 


■* 
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Input Card 6 

This card contains 4 curve fit coefficients for the upper bound 
of the Strouhal number vs. lambda/slgma function. They are as 
f ollovs : 


Word 1 

STUPK 

- +.25352^26+00 

Word 2 

STUFA 

- +.40487805+00 

Word 3 

STUPB 

- +.22229595+00 

Word 5 

STUPD 

- -.34329268-01 


Input Card 7 

This card contains 4 curve fit coefficients for the lower bound 
of the Strouhal number vs. lambda/slgma function. They are as 
follows: 


Word 1 

STL0K 

- +.1187C422+00 

Word 2 

STL0A 

- *.46569343+00 

Word 3 

STLOB 

• +.73139166-01 

Word 4 

STLOD 

- -.79927007-02 


Input Card 8 


This card contains 4 cu've fit coefficients for the critical curve 
of the Strouhal number vs. lambda/sigma function. They are as 
follows: 


Word 1 
Word 2 
Word 3 
Word 4 


PTCRITK - +.43502697+00 
STCRITA - - .61870504-01 
STCRITB - +.37269292-02 
STCRITD - +.40647482-02 


Input. Card 9 

This card contains 4 curve fit coefficients for the vortex force 
coefficient vs. lambda/slgma function. They are as follows: 


Word 1 CFK 
Word 2 CFA 
Word 3 CFB 
Word 4 CFD 


-.19458000+03 
+.25500000+02 
-.74460000+01 
-.39900Q0C K)0 


Input Cards 10 -through 15 

The curve fit coefficients for the Q-surface are read in at a rate o* 
four words per card, i.e. , QK (1), QA (1), 0B (1), QD (1) are punc' 
on Card 10; QK (2), QA (2), QB (2), QD (2) are on Card 11 of this. 
Reading continues per this format until JKAX »®ts of coefficients 
below have been read in. 


Card 10 QK(1) 

0 ( 1 ) 

V. ,1) 

QD(1) 

Card 11 QK(2) 

QA(2) 
QB(2) 
QD(2) 

Card 12 QK(3) 

QA(3) 

qbc: 

QD(3) 

Ca.:d 13 QK(4) 

QA(4) 
QB(4) 
QD(4) 

rd 14 QK(5) 

QA(5) 
QB(5) 
QD(5) 

Card 15 QK(6) 

QA(6) 
QB(6) 

v)D(6) 


- 4.0873881^+04 

- - 1.4052553E+02 

- 3.7419734E*r01 

- - 2.2574946E-03 

- J.3980471E+04 

- - 1.74986922+02 

« 3.8783556E+01 

- - 2.7034275E-03 

- 2 . 0081991E+04 

- - 1.4917770E+02 

- 4.5393842E+01 

- - 4.8689382E-03 

- 9.8799884E+03 

- - 1.2489887E+02 

- ..9950596E+01 

- - 6.8001116E-03 

- 7.8264710E+03 

- - 2. 068204 9E+OI 

- 4.3576094E+01 

- - 4.0612929E-03 

- 2.3506>69E+04 

- * 8.4432071E+02 

- 2.4773333E J -C1 

- - 1.4810690E-03 


Input C a »~d 16 


card contains exponent values (M) for material 1 


Word 

3 

XM 

(1,1) 

m - 

5.11 

(mean stress 

■i 

0 KSI) 

Word 

2 

XM 

(2,1) 

* - 

5.479 

( " 

ft 

m 

20 

" ) 

Word 

3 

XM 

(3,1) 

m - 

5.519 

( " 

ft 

m 

40 

" ) 

Word 

4 

XM 

(4,1 

- - 

5.645 

( " 

If 

m 

60 

" ) 

Word 

5 

XM 

(5,1) 

- - 

5.972 

( " 

ft 

m 

80 

" ) 


Input Card 17 

This card contains coefficient values (B) for material 1 


Word 

1 

B 

(1,1) 

- + .21410+16 

(mean stress 

m 

0 

KSI) 

Woi 1 

4- 

B 

(2,1) 

- + .72280+16 

( " 

It 

m 

20 

" ) 

Word 

3 

B 

(3,1) 

■ -f .44440+16 

( " 

fl 

*6 

40 

" ) 

Word 

4 

B 

(4,1) 

- + .24367+16 

( " 

ff 

* 

60 

" ) 

Word 

5 

B 

(5,1) 

- + .,£3200+16 

( " 

ft 

m 

80 

" ) 
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Input Card 18 


This card contains exponent values (M) for material 2 


Word 1 XM (1,2) 
Word 2 XM (2,2) 
Word 3 XM (3,2) 
Word 4 XM (4,2) 
Word 5 XM (5,2) 


13.003 

13.170 

16.008 

14.168 

5.345 


Input Card 19 


This card contains coefficient values (B) fox material 2 


Word 1 B (1,2) 
Word 2 B (2,2) 
Word 3 B (3,2) 
Word 4 B (4,2) 
Word 5 B (5,2) 


- + .67770+27 

- + .32560+28 

- + .29480+32 

- + .49910+27 

- + .89980+11 


Input Card 20 


This card contains exponent values (M) for material 3 


Word 1 XM 
Word 2 XM 
Word 3 XM 
Word 4 XM 
Word 5 XM 


(1,3) 

- - 2.447 

(2,3) 

- - 3.567 

(3,3) 

- - 4.387 

(4,3) 

* - 4.683 

(5,3) 

- - 6.124 


Input Card 21 

This card contains coefficient values (B) for material 3 


Word 1 B (1,3) 
Word 2 B (2,3) 
Word 3 B (3,3) 
Word 4 B (4,3) 
Word 5 B (5,*) 


- + .14360+10 

- + .13630+12 

- + .229004-13 

- + .12990+13 

- + .12510+13 


A. 6 Example Problem 


Listed below is an input data deck constructed in accordance with the 
instructions presented at the beginning of PROGRAM BELLOW. The notations 
that appear in columns 73 through 80 serve to identify the data group in each 
card. Following this listing is the corresponding computer output. The 
output is jrouped into three sections. The first group summarizes the perti- 
nent bellows input parameters. For this example, only the overall spring 
rate, KA, was inserted as data. The next group summarizes the fluid parameters. 
The next group contains the predicted longitudinal bellows performance. 

Bellows lock- in-range for a particular mode '*? vibration is defined by the 
upper and lower flow velocities. Stress indicator was calculated based on 


C - L 
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the critical flow velocity for each mode. Note, that for this particular 
bellows configuration, the lock- in-ranges for successive modes overlap, 
which indicates a more or less continuous spectrum of excitation velocities. 

Note also that all performance variables at the highest mode numbers are 
less than the corresponding quantities at previous mode numbers. Physically 
this behavior is accounted for t the fact that the apparent fluid mass is 
increasing at a faster rate than the dimensionless frequency numbers in Table A-I 
for t^is bellows. 
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B.l Flow Loop 

All liquid flow tests were conducted in a closed loop water flow 
tunnel shown schematically in Figure B-l. The bellows upstream piping was 
sized to the nominal bellows size, i.e. 3" PVC pipe was used during 3" 
bellows test and 6" PVC was used for 6" bellows. The flow loop can be 
pressurized to press*. res in excess of 100 pslg. 

Flow rate was accurately measured by a 4" turbine meter (Flow 
Technology SN - 64033 ) , and the loop's static pressure was determined by 
a calibrated bordon pressure gauge located one diameter upstream of the 
bellows . 


Fluid motion is generated by a Goulds propeller pump rated at 40 ft. 
head at 6000 GPM. The prime mover is a 75 hp variable speed hydraulic motor 
which provides a means to vary the loop flow velocity. Piping components are 
fabricated of carbon steel or PVC. A large antisurge reservoir (air over 
water) has been incorporated into the basic tunnel design. 

B.2 Bellows Instrumentation 

During a typical bellows flow test, three time dependent variables 
are normally recorded. These include the (1) volumetric flow rate, (2) the 
strain time history at various bellows locations, and (3) the displacement 
time history of selected bellows convolutes. 

The overall instrumentation setup is shown in Figure B-2 where it can 
be seen that three modes of recording data are possible. For quick look in- 
formation Polaroid pictures of the scope face may be obtained. As a second 
mode of operation, a high speed direct write galvonometer (CEC Model 5-124 ) 
is used to obtain high frequency hard copy bellows strain and displacement 
time histories; however, the majority of data (3rd mode of operation) was 
recorded in a form more useable for analysis, i.e. a dependent variable was 
plotted versus an independent variable on the x-y plotter while a test was in 
progress. 

Typical data collected in the form of two dimensional plots are presented 
in Figure 13. The vertical scale is proportional to either peak to peak strain 
amplitude or peak to peak displacement amplitude. Special circuitry, to be 
described subsequently, converted convolute peak to peak displacement motions 
to an equivalent D.C. analog voltage which was input to the y-axis of a model 
x-y recorder. Peak to peak strain signals (radial and circumferential) were 
procensed in a similar fashion. The horizontal axis is proportional to vol- 
umetric flow rate through the bellows. Since the primary flow measurement 
element was a turbine meter, its output frequency (directly proportional to the 
volume rate) was converted to a D.C. signal and then input to the recorder's 
x-axis . 


A typical instrumented bellows is shown in FigureB-3. Four strain gages 
were attached to the convolute crowns each test bellows. Two gages were placed 
on convolute number two, one responded to radial strains and the other responded 
to circumferential strain. Convolute number two was chosen ac a representative 
and convolute where peak strains occur (maximum relative displacement occurs in 
this region) but due to the end restraint. The middle convolute was gaged in 
the same manner as convolute number two. By observing the middle convolute 
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FIGURE B-l. WATER FLOW TUNNEL 
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Gage No . 


Gage Type 


Strain Direction Convolute No. 


EA-09-031ED-120 

EA-09-031ED-120 

EA-06-031DE-120 

EA-06-031DE-120 


Radial 2 
Radial 7 
Circumferential 7 
Circumferential 2 


FIGURE B-3 


STRAIN GAGE AND TAB LOCATION FOR 3" BELLOWS 
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response simultaneously with the second convolute, the mode number Is 
positively Identified and Insight Is gained with respect to the mode shape. 

All strain gages used were 1/32" long and each was selected to the 
base material of the bellows (321 stainless steel). Due to the small sire 
of the gage and its associated Installation difficulty, single arm active 
bridge circuits were employed. Figure B-4 shows a schematic of the signal 
conditioning circuit that converts gage resistance changes into a measurable 
voltage. The first stage amplifier (Analog Devices 610) is a high quality 
instrumentation amplifier operated in a differential voltage measurement mode. 
The second stage amplifier (Analog Devices 3140) provides offset voltage 
control and boosts the 610's output by a fixed gain of 10. 

Consolute displacement was obtained by measuring the displacement 
of a small metal tab that was epoxied to the crown of a convolute. A Bently 
probe (Model 316 ) was attached to a fixed structure above the test bellows. 
The tab couples with the transducer to produce an analog signal directly 
proportional to the displacement of the tab with respect to the transducer 
face; hence, the convolute absolute displacement was recorded. A sufficient 
number of tests were performed to insure that the virtually massless attached 
tab did not influence the vibration process. 




FIGURE B-4. STRAIN GAGE INSTRUMENTATION 
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BELLOWS GEOMETRIC AND MECHANICAL PROPERTIES DATA 



Bellows No 



15 

E 

3.3 

6.3 

3.0 

6.0 

3.6 

6.6 

. j 

.3 

.006 

.008 

2 

1 

13 

A 

.216 

.224 

.12 

.12 

1.8 

1.87 

93.5 

166.67 

42.37 

- 

3961 

> 
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E.l Two Phase Flow 


During the performance period of the bellows study, several 
special studies were conducted on an as needed basis. One parti- 
cularly noteworthy study conducted was a simplified analysis of 
the Shuttle LH 2 chilldown or recirculation system. Four possible 
operating conditions of the chilldown system were assumed and the 
analysis of the chosen "worst case," indicates low probability of 
a bellows failure due to a two phase flow phenomena. Results are 
presented below. 

E. 2 Case A - Pure Liquid Flow 

For this case the entire recirculation system was assumed 
to be flowing pure liquid hydrogen with the pump curve shown in 
Figure E.l defining the pressure head versus flow for each of the 
three pumps. Table E.l lists the assumed bellows geometry for 
this analysis. Table E.II lists the LH 2 properties and analysis 
results for the Case A pure liquid flow problem. As shown, because 
of the very low velocity and 1/2 pV 2 valve, the stress indicator 
valve is quite low and no bellows flow-induced vibration problem 
is anticipated. 

E. 3 Case B - Pure Gas Flow 


For this case we assume pure liquid flow through the pump 
followed by pure gaseous flow through the recirculation system. 

The reason for this assumption is to ensure the maximum possible 
driving head at the pump is available to "push" the gas through 
the lines. With gaseous flow through the pumps, a very low head 
would occur hence no means would exist to continue to introduce 
liquid into the system. 

It is assumed that sufficient heat is transferred into the 
liquid to cause complete boiling hence a pure gaseous flow through 
the recirculation lines. This is definitely a possiblitv at the 
first stage of chilldown. 



ps i 



0 0.5 1.0 1.5 2.0 

Flow Rate - lb/sec 


FIGURE E-l. RECIRCULATION PUMP PRESSURE 
FLOW CHARACTERISTIC 
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TABLE E-I Summary Of Bellows Data For Case A 


Bellows 

Geometry 

(Arrowhead Drawing 13619) 

• 

Material 

- 

ARMCO 21-6-9 


• 

Q 

• 

O 

- 

5.0 inches 


I.D. 

- 

4.0 inches 


N c 

= 

8 convo lutes 


N P 

= 

2 plys 


t 

r= 

0.008 inches per :-ly 


*c 

= 

2.0 inches convoluted length 


h 

= 

0.50 inches 


,\ 

= 

0.267 

• 

a 

= 

0.134 

Calculated Data 




k a 

s 

138.24 lb/inch overall spring rate 


M 

m 

= 

1.002 x 10” 4 lb-sec 2 /in 4 


f 

o 

= 

748 Hz, reference frequency 


f 

1 

= 

148.9 Hz, first mode frequency 


V 1 

= 

7.56 fps, first mode critical velocity 


f 

15 

= 

1488 Hz, highest longitudinal mode 
frequency 


V 15 


75.5 fps, highest longitudinal mode 
critical velocity 
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TABLE E -II Summary Of LH2 Properties And Analysis Results 
For Case A - Pure Liquid Flow 


Liquid Hydrogen Properties And Conditions 
LH 2 @ - 420°F, 16.1 psig 

p f = 0.002564 lbm/in 3 = 4.431 lbm/ft 3 

. (I) =4.2 lb/sec total flow, 3 pumps 


Calculated Data 


Volume flow = 


= 0.9479 ft /sec 


Volume flow 
Area 


0.9479 

0.08722 


= 10.87 fps 


1/2 p f V^ = 0.0565 psi 


CfQ = 8 (first mode) 

C F° h 2 2 

S.I. = (^f-) (£) (1/2 pV z ) 

N P 


= 882.8 psi 


Conclusions 


The Stress Indicator is so low no significant 
bellows response is possible. 
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Pure gaseous flow at the 4.2 lb/sec rate achieved for 
the pure liquid case is not possible because the flow loss 
would far exceed the available head at the pumps. There- 
fore, a downward adjustment in flow occurs until the loss 
matches the available pump head. The total flow from three 
pumps which satisfies this requirement is about 0.823 lb/sec; 
see Figure E-I. 

Based on this flow, the bellows of Table E-I has been 
analyzed and results are shown in Table E. III. As shown, the 
stress indicator is quite low and there is no possibility of 
acoustic resonance, hence the bellows is safe. 

E.4 Case C - Liquid Flow for Part of Line, Gaseous Flow 
for Rest 


For this case we assume pure liquid flow through the 
pumps and through a fraction of the total recirculation line 
length. The flow through the remainder of the recirculation 
line is assumed to be pure gaseous, '■'he transition from 
liquid to gas is assumed to occur suddenly at a single point 
in the line. 

As with Case B, the total pressure loss along the line 
is assumed equal to the head available from the pump. When 
the percentage of line with gaseous flow is large, we expect 
the mass flow to be smaller than the nominal 4.2 lb/sec value 
and the total head greater than the nominal 8.0 psi value. 

As a starting point we assume the presence of the gas will 
restrict the flow so that the pump is operating in the region 
of a 12 psi head value. Other assumptions are: 

. The liquid density is always p = 4.431 lb/ft^ 

c 

. The gas density is always = 0.0939 lb/ft^ 

. The pump head of 12 psi produces an average 

1/2 pV 2 of 0.0897 psi along the line 

. X is the percentage of the total li/.e over which 

the flow is pure liquid 

. (1-X) is the percentage of the total line length 

over which the flow is pure gaseous 

. There is sufficient heat transfer to convert the 

LH 2 to GH 2 at the point X 
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TABLE E-Ill. Summary Of GH? Properties And Analysis Results 
For Case B - Pure Gaseous Flow 


Gaseous Hydrogen Properties And Conditions 
GH 2 0 - 422°F, 19.0 psia 

p = 0.00005435 lbm/in 3 = 0.09392 lbm/ft 3 

y 

. d) » 0.823 lb/sec total flow, 3 pumps 

Calculated Data 

. Volume flow = — = 8.7t3 ft 3 /sec 

P 

g 

. V = 100.5 fps (could excite highest mode) 

1/2 pV 2 = 0.1023 psi 

. CfQ = 3.2 (highest mode) 

C F° h 2 2 

S.I. = {-£-) (£) (1/2 p v ) = 629.5 psi 

. Velocity required for acoustic resonance = 394 fps 

Conclusion 


Stress Indicator too low for problem. No acoustic 
resonance possible. Bellows safe. 
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We now have 

( I P L V l 2) X + ( I p g V g 2) (1 “ X) * °- 0897 P si 
From fluid continuity we find that 

<1 V’ * 47 -° P L V L 2) 

thus 

(1/2 p V 2 ) 

HTo — X + ( I p g V g } (1 " X) ■ 0.0897 psi 

or 

(1/2 p g V g 2 ) (1-0.0970 X) = 0.0897 psi 

From the above equation we find that a given value of X we 
have a unique vaiue of (1/2 pg V g 2) or V g . Figure E-2 shows 
a plot of Vg^ versus X from tne above equation. 

Note that as X increases toward a value of 1.0, the V g 
value also increases. For example, if there is liquid 
flow over 90% of the recirculation line, with the final 10% 
being gaseous flow, we can expect to have Vg = 272 fps from 
the gaseous flow over the final 10% of the line. 

Figure E-2 also shows the stress indicator values for the 
bellows defined in Table E-l. of course there must be a bellows 
located in the portion of the line over which the gaseous 
flow exists to experience this flow condition. 

From this analysis we can see that if the flow conditions 
assumed were to really exist then a bellows placed very near 
the end of the recirculation line could be subject to rather 
high stresses. Also we are getting into gaseous velocity 
ranges where acoustic resonances might be possible. Table E-IV 
summarizes the results of this analysis. 

The question remaining then is: Can such a flow condi- 

tion occur? The answer to this question depends on the re- 
sults of a heat transfer analysis to find out if sufficient 
heat can be introduced into the fluid to produce the required 
phase change from liquid to gas. 


6,000 



FIGURE E-2. RESULTS OF CASE C ANALYSIS 
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TABLE E- IV Case C - Portion Of Line Pure Liquid Flow 
And Portion Pure Gaseous Flow 


Hydrogen Properties And Conditions 
LH 2 § - 420°F, 16.1 psig 

P L « 4.431 lbm/ft 3 

GH 2 = -422°F , 19.0 psia 

p g = 0.09392 lbm/ft 3 

. ' Fluid mass flow variable 

Analysis 

. Pure liquid flow over X percent of line length 

. Pure gaseous flow over (1-X) percent of length 

2 

. Pump head always 12 psi, average line 1/2 pV = 
0.0897 psi. 

. X and flow head related by 

(1/2 p„ V 2 ) (1-0.979X) = 0.0897 psi 

g g r 

. Solution to above given in Figure 2 

. For example, if X =.90 or 90%, then 

1/2 p„ V 2 = 0.754 psi 

g g 

V = 272.6 fps 

9 

S.I. = 4713 psi (SAFE) 

. Acoustic resonance occurs ® V * 394 fps 

9 

Conclusions 


Only a bellows located at end of line would be in 
possible damage if the postulated flow condition 
can actually occur. 
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This question will be answered in the next section; however, 
if no bellows exists over the final 10% of the line length 
then no problem exists. 

E. 5 Case D - Slug Flo.: 


For this case we assume that a pocket of gas has formed 
in the recirculation line and is growing because of further 
boiling of LH 2 . This gas product growth pushes the LH 2 in 
front of it out of the line; hence, we need to determine if the 
liquid and/or gas velocities can become high enough to create 
a bellows problem. 

Figure E-3 shows a schematic diagram of the physical orob- 
lem for Case D. As shown, we assume a gas pocket of length 
Y which is growing because of boiling caused by heat trans- 
fer through the tube from the surroundings. The rear bound- 
ary of the gas pocket is assumed moving at a velocity Vx, 
while the front boundary is assumed moving at a velocity V 2 . 

The difference in velocity of the two boundaries relates to 
the volume growth of the gas pocket because of boiling. 

The first problem to be solved is a determination of the 
boiling volume growth of the gas pocket. Table E-V summarizes 
an analysis to solve this problem. We assume a gas pocket 
of length Y is being f o riuC d by boiling from heat transferred 
through the tube wal.l. It has been determined that the 
boiling transfer coefficient on the inside of the tube is so 
very high relative to the external heat transfer coefficient 
that the tube wall can be assumed at the same temperature 
as the LH 2 . Therefore the boiling rate is limited or deter- 
mined by the heat transfer from the ambient surroundings to 
the tube wall. 

On this basis the analysis shows that the maximum weight 
rate of LH 2 boiled into GH 2 will be 


u> = 5.29 lbm/hr 
max 

per foot of tube over which boiling is assumed to occur. 
From this rate of boiling the volume rate of growth of the 
gas pocket has been calculated to be 

Q , * 0.9199 ft^/min per foot length, 
vol 
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FIGURE E-3, Schematic of Cese D Problem 



E-12 


TABLE E-v Summary Of Heat Transfer Through 
Recirculation Line Walls 



Q = heat transfer through wall to induce boiling 
of LH 2 

Y = length of line over which boiling assumed 
occurring 

LH 2 assumed @ - 420°F, 16.1 psig 

transfer iimited by convection to. tube on O.D. - 
wall assumed at temperature virtually equal to 


above 


Heat 

tube 

lh 2 

From 


Q = h 0 A (T a - T t ) 

h Q = convection neat transfer coefficient 
assumed equal to 2.0 Btu/hr-ft2 °F 


A = ttD^Y = area of tube O.D. for length Y 
Per foot of tube we have where T = 70°F and T_ = -420 

cl 1 

Q = 1026 Btu/hr per foot of tube 

If t.ie heat of vaporization of LH~ is assumed at 
194 Bt' /lbm then the weight rate of fluid boiled is 


01 


1026 Btu/tr , n_ ^ 

‘ l§4 ~ B tu /ib " m = U29 lbm/hr P er £oot 


of tube 

From above the relative boundary velocities of the 
gas pocket has been calculated at 


V 


2 


0.1758-Y fps 
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Finally this volume growth rate permits calculation of the 
differential gas pocket boundary velocities as 


V 2 - V 1 * 0.1758 Y-fps 

From the above it is clear that boiling over very long 
lengths of line would be required to cause significant 
increases in the advancing liquid-gas boundary. For example, 
we might make some probably impossible assumptions to show 
that there is no real problem from bellows flow excitation 
for the Case D situation. Let's assume: 

. The liquid weight flow at the rear boundary is 

<I>£ = 4.2 lb/sec. 

. The rear boundary advances at a rate correspond- 

ing to the above or, V. = 10.87 fps (see Table 
E-II ) . 1 

. Boiling occurs ov<. r a 50 foot length of recircu- 

lation line. The line may or may not be this 
long . 

Based on the above, we have: 

V 2 = 10.87 + 0.1758 x 50 - 19.66 fps 

The liquid in front of the gas boundary is therefore being 
"pushed" along at a velocity of 19.66 fps. The stress 
indicator for this particular case would be (CfQ = 2.67, 

3rd mode) 

S . I . = 963.7 psi 

which is clearly too low to cause any problem. 

Discussion and Conclusions 

Figure E-4 shows a realistics but simplified schematic of 
the LH 2 feed and recirculation systems. During chilldown the 
prevalves are closed, the recirculation pumps are operative and 
the recirculation valves are open. From our analysis so far we 
anticipate the following chain of events. 

(1) LH 2 will start to flow into the feed system from 
the recirculation pumps at a rate greater than 
the nominal 4.2 lb/sec since the system is empty. 

(2^ Massive boil off will initially occur as the feed 
lines and pumps begin to cool down. 


I 
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(3) The initial boil off wi. 1 1 raise the gas pres- 
sure in the feed lins and pump areas, but as 
the pressure increases the LH 2 flow from the 
recirculation pumps will slow down or shvt off 
as the maximum pump head pressure is achie ed. 


(4) The initial flow through the recirculation . ine 
will be pure gaseous under conditions outlined 
in Case B. 

(5) As the feed system and pumps begin to cool 
down, LH 2 will enter the recirculation lines. 

We can expect a condition of slug flow where 
we have alternate pockets of gas and liquid. 

As shown in Case D, there is not sufficient 
heat transfer into the recirculation lines to 
create a high velocity condition from local 
boiling. 

(6) When the system is chilled down to the required 
extent, pure liquid flow will occur and Case A 
analysis covers this situation. 

The Case C analysis is, we feel, unrealistic since, as 
shown in the Case D analysis, we cannot expect sufficient 
heat transfer through the recirculation lines to achieve boil- 
ing at a rate necessary to create a high velocity problem. 

At this time we feel there is little chance of a bellows 
related problem in the feed and recirculation system because 
of two phase flow problems. 



